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THE  FORCES  AND  MOMENTS  ON  A  MISSILE  DUE  TO  A 
FLUID  JET  EJECTING  LATERALLY 


Volume  II.  Computer  Program  and  Operating  Instructions  (SIDJET) 


1.  INTRODUCTION 

Under  Contract  No.  N00024-79-C-5381  with  Naval  Sea  Systems  Command, 
Nielsen  Engineering  &  Research,  Inc.  (NEAR)  has  developed  engineering 
aerodynamic  prediction  methods  for  calculating  lateral  jet  induced 
forces  and  moments  on  f inless  missiles,  missiles  with  monoplane  wings 
and  on  missiles  with  cruciform  fin  configurations.  Those  methods 
concern  the  problem  of  accounting  for  the  forces  and  moments  induced  by 
interference  between  jets  issuing  from  the  side  of  a  missile  and  the 
external  flow  around  the  missile.  For  bodies  with  no  fins,  the  method 
of  least  squares  was  used  to  correlate  the  available  data  to  obtain 
simple  formulas  for  the  jet  induced  force  and  moment.  For  cruciform 
finned  bodies,  an  extensive  set  of  data  from  wind  tunnel  tests  was  used 
to  develop  a  tabular  data  base  from  which  the  forces  and  moments  can  be 
calculated.  For  bodies  with  monoplane  wings,  a  combination  of  the 
body-alone  formulas  and  the  cruciform  data  is  used. 

This  report  presents  instructions  for  using  the  computer  program 
developed  to  implement  the  prediction  methods.  The  program  is  constructed 
as  a  primary  subroutine  (SUBROUTINE  SIDJET)  and  several  auxiliary  sub¬ 
routines.  The  calling  program  can  be  simply  a  program  for  reading  the 
required  input  data  and  printing  the  output  quantities  or  a  more  sophis¬ 
ticated  program  of  any  kind. 

The  remainder  of  this  report  describes  the  equations  and  procedures 
programmed,  the  input  requirements  of  the  program,  and  the  output 
quantities  that  are  calculated  by  the  program.  Sample  cases  are 
presented  to  aid  in  checking  out  the  program. 


1 


2.  EQUATIONS  AND  PROCEDURES  PROGRAMMED 


2.1  Forces  and  Moments  on  a  Body  of  Revolution  With  a  Radial  Jet 

The  general  problem  is  of  jets  issuing  from  a  missile  body  and  inter¬ 
fering  with  the  body  pressure  distribution.  During  the  outward  movement 
of  the  jet,  it  turns  downstream  as  a  result  of  interacting  with  the  free 
stream.  During  its  downstream  passage  the  jet  entrains  free-stream  air. 

If  the  jet  exhausts  into  a  vacuum,  the  reaction  force  on  the  missile 
is  dependent  only  on  the  jet  momentum.  However,  because  there  is  airflow 
around  the  missile,  interference  between  the  jet  and  the  external  airflow 
can  cause  an  additional  external  reaction  which  can  either  amplify  or 
reduce  the  net  force  on  the  missile  and  create  an  additional  moment  due 
to  the  displacement  from  the  jet  centerline  of  the  center  of  pressure  of 
the  flow  field. 

The  amplification  of  the  jet  reaction  force  is  defined  as 


where 


PT  +  Fi 
^  =^F-i 


(1) 


F  =  the  jet  thrust  without  external  flow 
T 

=  the  additional  reaction  due  to  the  jet/ external  flow 
interaction 

The  moment  produced  by  the  jet/external  flow  interaction  is 
defined  in  terms  of 


where 


F .  (x  .  -  x  .  ) 
l  j _ i_ 

ftd 


(2) 


the  jet  axial  station  measured  positive  aft 
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=  the  axial  location  of  the  center  of  pressure  of  the 
interaction  force 

D  =  the  body  (reference  or  maximum)  diameter 

The  expressions  used  to  calculate  and  (derived  from  least  squares 

fits  of  experimental  data)  are  ^ 


=  0.6118  +  0.1358(1  -  0.485  /T75)//C~ 
+  0.0946Moo  +  0.004317/U/D) 


and 

=  0.5582  -  0.1884//CZ  -  1.9659/U/D)  (4) 

Mj  T 

where 


£  =  axial  distance  of  jet  from  base  of  missile 

CT  =  jet  thrust  coefficient,  F^/q^S 
M  =  free  stream  Mach  number 

oo 

q^  =  free  stream  dynamic  pressure 

S  =  reference  (maximum)  cross  sectional  area  of  body 

The  forces  and  moments  corresponding  to  equations  (3)  and  (4)  for  a 
particular  configuration  are  calculated  using  the  right  handed  coordinate 
system  shown  in  the  sketch. 
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Jet 


The  normal  force  and  moment  coefficients  in  the  plane  of  the  jet  are: 


CN  =  -CTKN 


(5) 


CM .  CTKM . 
J  D 


(6) 


and  the  components  of  the  force  and  moment  in  the  coordinates  shown  are: 


C  = 


C„  = 


m . 
3 


n . 
3 


C7cos  <i> . 
N  Yj 


C  sin  (j>  . 
N  Yj 


CM.C°S 

CM.Sin 


C„  =  0 


(7) 

(8) 

(9) 

(10) 
(11) 


2 . 2  Forces  and  Moments  on  a  Cruciform  Wing-Body  Combination  with  a 
Radial  Jet 

Addition  of  wings  and  fins  to  the  bodies  of  revolution  discussed 
previously  introduces  major  new  complications  in  estimating  the  forces 
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and  moments  resulting  from  the  jet  issuing  from  the  body.  Most  obvious, 
the  fins  permit  development  of  rolling  moments  when  acted  on  by  the  pres¬ 
sure  field  and  trailing  vortices  of  the  jet.  The  fins  also  allow  major 
increments  in  the  normal  and  side  forces  and  their  accompanying  trans¬ 
verse  moments . 

The  extreme  complexity  of  the  jet-body-fin  interaction  problem 
precluded  the  development  of  simple  formulas  for  the  forces  and  moments . 
Consequently,  recourse  was  made  to  a  data  base  in  which  the  five  force 
and  moment  coefficients,  ,  G^,  Cmj ,  and  are  tabulated  as 

functions  of  the  missile  angle  of  attack,  a,  and  roll  angle  (f>m  (see 
sketch),  the  jet  axial  position,  x.  and  the  jet  azimuthal  position,  <J>  .  . 


The  roll  angle,  <f>m,  is  the  angle  between  the  plane  of  the  wing  (number  1 
in  the  sketch)  which  is  adjacent  to  and  clockwise  from  the  jet  (when 
viewed  from  the  rear)  and  the  plane  containing  the  missile  axis  and  the 
flight  velocity  vector.  The  angle  $  is  positive  clockwise.  The  jet 
azimuthal  angle,  <Jk  ,  is  the  angle  between  the  jet  axis  and  wing  number  1. 
The  angle  cj^  varies  from  0°  when  the  jet  blows  along  wing  number  1  to 
-90°  when  the  jet  blows  along  wing  number  4.  The  angle  of  attack  varies 
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from  -10°  to  +10°  in  the  data  base.  The  jet  axial  position  is  tabulated 
in  terms  of  the  position  relative  to  the  fin  (wing)  root  leading  edge, 

(Xj  -  XLE^LF'  wlier®  is  tlie  (win9)  root  chord  and  xLE  is  the  axial 

coordinate  of  the  fin  (wing)  root  leading  edge.  The  force  and  moment 
coefficients  are  normalized  by  the  jet  thrust  coefficient.  A 
proportional  relationship  between  the  interaction  forces  and  moments  and 
the  jet  thrust  is  used,  based  on  experimental  measurements  at  three 
different  thrust  levels.  The  free  stream  Mach  number  of  the  experimental 
tests  was  1.6  for  most  of  the  data.  A  few  tests  at  Mach  numbers  of  1.8 
and  2.1  indicated  that  a  reasonable  first  approximation  is  that  Mach 
number  effects  over  this  range  are  negligible.  Therefore,  variations 
with  Mach  number  are  not  included  in  the  tabulated  data  base.  It  is 
felt  that  the  data  base  can  be  used  in  the  low  supersonic  Mach  range 
with  reasonable  results,  but  it  has  not  been  validated  for  subsonic 
speeds  because  of  the  lack  of  data.  Likewise,  variations  with  fin  (wing) 
span  are  not  included  because  such  effects  were  not  separable  from  the 
effect  of  the  chordwise  position  of  the  jet  for  the  swept-leading-edge 
fin  used  in  the  wind  tunnel  tests.  The  values  of  the  independent 
variables  used  in  the  tabulated  data  are  as  follows: 


a  :  0°,  ±2°,  ±4°,  ±10° 


0.19,  0.54,  0.89 


:  -22.5°,  -45°,  -67.5° 


cf>m  :  0,  22.5°,  45°,  67.5°,  90°,  135°,  180°,  202.5° 

225°,  247.5°,  270°,  315°,  360° 


The  tabulated  data  are  interpolated  in  the  following  manner:  First, 
the  input  values  of  a,  (f>m,  <J> ^  and  x^  are  compared  with  the  values  corres¬ 
ponding  to  the  tabulated  data  to  obtain  bracketing  values  for  each 
quantity.  Those  sixteen  quantities  are  reduced  to  eight  by  interpolating 
to  the  input  value  of  a.  The  eight  interpolated  values  are  reduced  to 
four  by  interpolating  in  $  .  The  four  new  values  are  reduced  to  two  by 
interpolating  in  <j> ^  ,  and  the  final  interpolation  in  x^  produces  the  result. 
For  values  of  the  input  parameters  outside  the  range  of  the  tables,  the 
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results  should  be  used  with  caution.  Experimental  measurements  were 
available  for  angles  of  attack  up  to  approximately  +  12°,  but  the 
variation  of  the  interaction  forces  and  moments  for  higher  angles  of 
attack  is  not  known.  Likewise,  experimental  data  were  obtained  for 
values  of  <(>j  of  -67.5°  and  -45°,  with  the  tabulated  data  at  <j>j  =  -22.5° 
obtained  by  symmetry  relations.  The  effect  of  a  jet  nearer  a  fin  than 
22.5°  is  not  expected  to  be  drastically  different  than  at  22.5°.  In 
order  to  provide  the  most  reasonable  results  for  analyzing  the 
trajectory  of  a  vehicle,  the  boundary  values  of  the  tabulated  data 
are  used  for  input  parameters  outside  the  range  of  the  tables  rather 
than  an  extrapolation,  because  calculations  showed  extrapolated 
values  to  be  quite  unreasonable  in  some  cases. 

2 . 3  Forces  and  Moments  on  a  Monoplane  Wing-Body  Combination  With 
a  Radial  Jet 

The  calculative  procedure  for  a  monoplane  winged  configuration  uses  a 
combination  of  the  procedures  for  cruciform  configurations  and  body- 
alone  configurations.  To  illustrate  the  procedure,  the  sketch  below 
shows  a  missile  viewed  from  the  rear  at  a  small  roll  angle.  This  is 
essentially  the  same  configuration  as  was  shown  for  cruciform 
configurations  except  that  fins  1  and  3  are  assumed  to  be  missing. 
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The  induced  normal  force  and  pitching  moment,  Cz  and  Cm .  on  the 
cruciform  configuration,  are  produced  primarily  by  the  jet  thrust  and 
the  induced  force  on  the  body  and  fin  number  4.  For  the  monoplane  case, 
it  is  assumed  that  the  removal  of  fin  number  1  has  a  negligible  effect 
on  the  manner  in  which  the  jet  affects  the  flow  over  fin  number  4. 
Likewise,  the  removal  of  fin  3  has  a  negligible  effect  on  the  manner  in 
which  the  jet  influences  fin  number  2,  particularly  since  the  two  fins 
are  on  the  side  away  from  the  jet.  Thus,  the  Cz  and  Cm .  for  a 
monoplane  configuration  are  determined  directly  from  the  cruciform  data 
base . 

In  the  absence  of  the  vertical  wing  the  induced  side  forces  can  only 
be  produced  by  the  jet  and  its  interaction  with  the  body.  It  is  assumed 
that  the  influence  of  the  horizontal  wing  on  the  jet  is  negligible. 
Therefore,  the  incremental  side  force  and  moment  are  taken  to  be  those 
for  the  case  of  a  jet  from  a  body  without  wings.  Thus, 

Cy  =  Cy  for  a  body  without  wings 

for  a  body  without  wings 

The  recommended  relationship  for  the  induced  rolling  moment  is 

Cff  =  +0.6  C„  sin  d>  . 

£  N 

in  which  the  term  sin  <J>  .  has  been  inserted  in  order  to  obtain  zero  roll 

.  3 

moment  for  a  jet  flowing  in  the  plane  of  symmetry. 

2.4  Forces  and  Moments  on  a  Body  with  Both  Wings  and  Tails 

No  provision  is  made  in  the  program  for  calculation  of  forces  and 
moments  on  a  body  with  both  wings  and  tails.  Such  calculations  can  be 
made,  however,  by  simply  providing  the  required  information  for  the 
set  of  surfaces  that  are  most  likely  to  interact  with  the  jet.  For 
example,  if  the  jet  is  located  ahead  of  the  wing  trailing  edge,  the 
wings  themselves  would  probably  be  the  most  likely  surfaces  to  be 
defined  in  the  input  to  the  program.  On  the  other  hand,  if  the  jet  is 
located  aft  of  the  wing  trailing  edge  but  ahead  of  the  tail  trailing 
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edge,  the  tail  surfaces  should  be  input.  Some  judgment  must  be  exercised 
by  the  user  in  determining  the  appropriate  input  quantities  depending  on 
the  relative  sizes  of  the  wing  and  tail  surfaces  and  their  proximity  to 
each  other  as  well  as  to  the  jet. 


3.  PROGRAM  OPERATING  INSTRUCTIONS 


3 . 1  Input  Data  Required 

Communication  with  subroutine  SIDJET  is  through  a  formal  parameter 
list.  Programming  of  specific  input/output  procedures  is  left  for  the 
user  to  provide.  The  input  quantities  required  by  the  program  are  as 
follows : 


Variable  Name 


Description 


NCONFG 


XJIN 

PHIJIN 

PHIMIN 

ALPHIN 

L 

D 


An  index  for  identifying  the  type  of  configuration  under 
consideration 
0  -  body  alone 
2  -  monoplane 
4  -  cruciform 

Any  other  value  input  will  cause  the  printing  of  the 
message  "configuration  index  in  error" 

Distance  of  jet  from  nose 

Azimuthal  location  of  jet,  4>  j  /  degrees 

Roll  angle  of  missile,  <|>  ,  degrees 

Angle  of  attack,  a,  degrees 

Total  length  of  missile 

Maximum  or  reference  diameter  of  missile 
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CT 


MINF 

XLE 


Jet  thrust  coefficient,  =  jet  thrust  with  no 
external  f  low/ ( q^TrD^/4 )  .  This  can  be  approximated  by 


where  pa  =  ambient  static  pressure 
pt  =  jet  stagnation  pressure 
Yj  =  jet  ratio  of  specific  heats 
Flight  Mach  number,  (.6-4 .0) 

Distance  from  nose  to  fin  root  leading  edge 


LF 


Length  of  fin  root  chord 


OPFLAG  Initialization  index;  should  be  set  to  zero  in  the  main 

program  before  the  first  call  to  subroutine  SIDJET.  In 
normal  operation,  OPFLAG  is  then  set  to  999  in  subroutine 
TABULA  after  initialization  of  certain  arrays.  Subsequent 
calls  to  SIDJET  with  OPFLAG  unchanged  cause  the  initializa¬ 
tion  to  be  passed  over. 


3.2  Special  Instructions  for  Preparing  Input  Data 

3.2.1  Units . -  All  quantities  with  length  dimensions  can  be  input  in 
any  units  including  nondimensional .  The  main  rule  is  that  all  units 
must  be  consistent.  If  one  of  the  input  lengths  is  chosen  as  a  reference 
quantity  to  nondimensionalize  the  other  lengths,  the  reference  quantity 
should  be  input  as  1.0.  All  lengths  are  subsequently  nondimensionalized 
in  SIDJET  using  either  the  body  diameter  or  the  fin  length  independent 
of  any  previous  nondimensionalization  of  the  input  data.  All  angles  are 
input  in  degrees. 
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3.2.2  Angles  . -  The  two  angles,  PHI JIN  and  PHIMIN  require  special 
consideration  of  the  type  of  configuration  being  treated.  Different 
coordinate  systems  are  used,  depending  on  the  configurations.  The 
finless  case  (body  alone)  uses  the  crossflow  wind  as  reference.  The 
monoplane  uses  the  "top",  since  such  configurations  are  bank-to-turn , 
aircraft-like  vehicles.  The  cruciform  configuration  has  no  preferred 
bank  angle,  so  the  fin  adjacent  to  the  jet  is  used  as  a  reference. 


Finless  Configuration 


Measure 


clockwise  from  a  leeward  reference 


plane  (looking  forward  from  the  tail  of  the 
missile) .  The  allowable  range  of  ^  is 
0  _<  (fj  <  360°.  The  roll  angle,  ^  is  not  used 
for  this  configuration.  The  forces  and  moments 
are  resolved  into  components  in  the  right-handed 
system  shown. 


z 


Monoplane  Configuration 


Designate  the  top  side  symmetry  plane  of  the  vehicle  as  reference. 
Measure  <j>  .  clockwise  from  that  plane  to  the  jet  axis.  Measure  <j> 
clockwise  from  the  leeward  (for  positive  angle  of  attack)  plane  con¬ 
taining  the  wind  vector  and  the  body  axis  to  the  reference  plane. 

Both  <j>  .  and  d>  can  be  between  0°  and  360°.  The  forces  and  moments  are 
D  m 

then  resolved  into  the  y-z  coordinate  system  shown 


Top  side 
reference 
plane 
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Cruciform  Configuration 


The  cruciform  fin  data  tables  apply 
to  a  jet  located  between  two  particular 
fins.  Thus,  for  a  given  jet  position, 
the  fin  nearest  the  jet  clockwise  from  (3) 
the  jet  should  be  designated  as  fin 
number  1 .  The  other  fins  are  numbered 
clockwise  from  fin  number  1.  The  angle  <j> 
is  measured  from  fin  number  1  to  the  jet 


( -90c 


<  0°) .  The  angle 


IS 


—  J  —  Tm 

measured  from  the  wind  vector  to  fin 

number  3  as  shown  (for  positive  angle 

of  attack) .  The  forces  and  moments  are 

determined  in  the  y-z  coordinate  system 

shown,  that  is ,  in  general  z  is 

positive  along  fin  number  1  and  y  is 

positive  along  fin  number  2.  All  moments 

are  referred  to  the  jet  location. 


3 . 3  Output 

The  output  from  program  SIDJET  consists  of  the  force  and  moment 

coefficients  as  follows: 

CZ  component  of  normal  force  in  the  z  direction  (see 

sketches  in  previous  section  for  definition  of  coordi¬ 
nate  axes  for  the  various  allowable  configurations) 

CY  component  of  normal  force  in  the  y  direction 

CMJ  component  of  moment  about  the  y  axis  at  the  jet  . 

axial  position  -  positive  nose  up 

CNJ  component  of  moment  about  the  z  axis  at  the  jet 

axial  position  -  positive  nose  right 

CL  rolling  moment  about  the  missile  axis  -  positive 

right  wing  down  or  clockwise  looking  forward 
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4. 


SAMPLE  CASES 


4.1  Wingless  Body 

Input  quantities  required  for  a  sample  calculation  on  a  wingless 
body  are  as  follows: 

NCONFG  =  0 
XJIN  =  4.364 
PHIJIN  =  20° 

PHIMIN  =  0° 

ALPHIN  =  0° 

L  =  9.0 

D  =1.0 

CT  =  .4068 

MINF  =1.5 

XLE  =0.0 

LF  =0.0 

OPFLAG  =  0 


The  output  quantities  will  then  be 

Cz  =  -.2849 
Cy  =  -.1037 
CMJ  =  -.0616 
CNJ  =  -.0224 
CL  -  0 


4.2  Cruciform  Wing-Body  Configuration 

Input  quantities  for  a  sample  calculation  on  a  cruciform  wing-body 
combination  are  as  follows: 

NCONFG  =  4 
XJIN  =20.0 
PHIJIN  =  -40.0 
PHIMIN  =30.0 
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ALPHIN 

= 

o 

i 

L 

= 

35.63 

D 

= 

2.52 

CT 

= 

1.0 

MINF 

= 

1.6 

XLE 

= 

15.54 

LF 

= 

17.41 

OPFLAG 

= 

0 

The  output  quantities  will  be 

CZ  =  -.0621 
CY  =  -.0282 
CMJ  =  -2.5317 
CNJ  =  1.9349 
CL  =  .0092 


4.3  Monoplane  Wing-Body  Configuration 

Input  quantities  for  a  sample  calculation  on  a  monoplane  wing-body 
combination  are  as  follows: 


NCONFG 

= 

2 

XJIN 

= 

20.0 

PHIJIN 

= 

320.0 

PHIMIN 

= 

30.0 

ALPHIN 

= 

-4.0 

L 

= 

35.63 

D 

= 

2.52 

CT 

= 

1.0 

MINF 

= 

1.6 

XLE 

= 

15.54 

LF 

= 

17.41 

OPFLAG 

= 

0 
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The  output  quantities  will  be 

CZ  =  -.0621 
CY  =  .4748 
CMJ  =  -2.5317 
CNJ  =  .0271 
CL  =  .0240 


5.  COMPUTER  PROGRAM  LISTING 


SUBROUTINE  SIDJET  (NCONFG.XJIN.PHI JIN.PHIMIN, ALPHIN.L.D. 

*  CT.MINF.XLE.LF.OPFLAG.CN.CY.CMJ.CNJ.CL) 

REAL  L.LD.MINF.LF 
INTEGER  OPFLAG 

THIS  SUBROUTINE  IS  DESIGNED  TO  CALCULATE  NORMAL  Y  AND  Z  FORCES 
AND  PITCHING.  ROLLING  AND  YAWING  MOMENTS  OF  A  MISSILE  GIVEN: 


INPUT  VALUE 

VARIABLE  NAME 

RANGE 

CONFIGURATION 

NCONFG 

0,2,4 

JET  LOCATION 

XJIN 

JET  AZIMUTHAL  POS. 

PHIJIN 

-90.0  TO  0.0 

MISSILE  ROLL  ANGLE 

PHIMIN 

0  TO  360 

ANGLE  OF  ATTACK 

ALPHIN 

-10  TO  10 

THRUST  COEFFICIENT 

CT 

MISSILE  LENGTH 

L 

MISSILE  DIAM.  (REF.) 

D 

AXIAL  LOCATION  OF 

WING  ROOT  LEADING 
EDGE 

XLE 

WING  ROOT  CHORD 

LF 

MACH  NUMBER 

MINE 

TABLE  INITIALIZATION 

OPFLAG 

FLAG  SHOULD  BE  C 

SET  TO  ZERO  IN  MAIN 
AT  START  OF  PROGRAM 


THIS  IS  DONE  EITHER  THROUGH  CALCULATION  OF  THESE  FORCES  AND  MOMENTS 
OR  THROUGH  INTERPOLATION  OF  TABLES  OF  EXPERIMENTAL  AND  CALCULATED 
DATA.  DEPENDING  ON  THE  CONFIGURATION  OF  THE  BOOY  IN  QUESTION.  I.E.* 
WINGLESS,  MONOPLANE,  OR  CRUCIFORM  C 
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c  C 
C  THIS  SUBROUTINE  RETURNS  THE  FOLLOWING  VALUES:  C 
C  c 
C  OUTPUT  VALUE  OUTPUT  NAME  C 

c  — -  -  c 

C  NORMAL  FORCE  COEFFICIENT  CN  C 
C  SIOE  FORCE  COEFFICIENT  CY  C 
C  PITCHING  MOMENT  ABOUT  CMJ  C 
C  JET  AXIAL  POSITION  C 
C  YAWING  MOMENT  ABOUT  CNJ  C 
C  JET  AXIAL  POSITION  C 
C  ROLLING  MOMENT  CL  C 
C  C 
C  C 


COMMON  /MASTBL/  DCNX 1(13*3*7)* DCNX2 (13*3*7) , DCNX3 (13.3.7) » 

*  DCYXi ( 13,3,7) »DCYX2(13.3.7) .DCYX3 ( 13*  3*7) • 

*  DCMX1 (13*3*7) .DCMX2( 13.3.7) .DCMX3 ( 13.3.7)  ♦ 

*  DCHX1 (13*3*7) .DCHX2 (13*3*7) .DCHX3 < 13.3.7) * 

*  0CLX1 (13.3.7) ,DCLX2( 13.3.7) .DCLX3 ( 13.3.7) 
I80UND=0 

JUMPMP=0 

C 

C  DECIDE  ON  THE  BASIS  OF  CONFIGURATION  NUMBERS  WHETHER  THE 
C  MISSILE  IS  OF  THE  BODY-ALONE.  CRUCIFORM  OR  MONOPLANE  TYPE 

IF  BODY  ALONE  CONFIGURATION,  USE  CALCULATED  VALUES  OF  FORCE 
AND  MOMENT 

IF  (NCONFG.EO.O)  GO  TO  100 

IF  WINGED*  BUT  JET  TRAILS  ALL  WINGS  AND/OR  TAILS.  USE  BODY-ALONE 

IF  (NCONFG.EQ.2. AND. XJIN.GT .XLE+LF)  60  TO  100 
IF  (NCONFG. EQ. 4. AND. XJIN.GT. XLE+LF)  GO  TO  100 

IF  CRUCIFORM  CONFIGURATION.  VALUES  OF  FORCE  AND  MOMENT  ARE 
INTERPOLATED  FROM  TABLES  OF  EXPERIMENTAL  AND  CALCULATED  DATA 

IF  (NCONFG. EQ. 4. AND. XJIN.LE.XLE*LF>  GO  TO  200 

IF  MISSILE  IS  MONOPLANE,  AND  JET  IS  NOT  AFT  OF  WING  SURFACES. 

A  COMBINATION  OF  BODY-ALONE  AND  CRUCIFORM  CALCULATIONS  IS  DONE 

IF  (NCONFG. EQ. 2. AND. XJIN.LE.XLE*LF)  GO  TO  300 

IF  THE  CONFIGURATION  IS  NONE  OF  THESE,  EXIT  THE  SUBROUTINE 

60  TO  400 

CALCULATE  NORMAL  FORCES  AND  MOMENTS  FOR  BODY  ALONE  CONFIGURATION 
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100  LD=(L-XJIN)/0 

IF  (LD.GT.6.0)  LD  =  6.0 
IF  (CT.LT.0.0.OR.LD.LE.0.0)  GO  TO  120 

CNT=0.611845*  (0.135753/SQRKCT)  > * ( 1 . 0- ( 0 ,484984*SQRT (LD)  )  ) 
*  ♦ (0.094633*MINF) ♦ (0.004317/LD) 

CNT=-CNT 

CMJ1=. 558196- ( . 188399/SQRT (CT ) ) - ( 1 .96586/LD) 

GO  TO  121 

120  CNT  =  0.0 
CMJ1  =  0.0 

121  PHIRAD=PHIJIN*. 0174533 

IF  (JUMPMP.NE.O)  GO  TO  125 
CN=CNT*COS(PHIRAD)*CT 
CY=CNT*SIN(PHIRAO)*CT 
125  IF  (CT.EQ.0.0)  GO  tO  130 

XI  =  XJIN-D*(CMJl/(  l-ABS(CNT)  )) 

IF  (XI.GT.L)  CMJ1  =  -  <  ABS (CNT)  -1 ) * ( (XJIN-L) / ( 2*D) ) 

CMJ1=CT*CMJ1 

IF  (JUMPMP.NE.O)  GO  TO  355 
CM J  =  CM J 1  *  COS ( PH IRAQ ) 

60  TO  132 
130  CM J=0 . 0 

132  CNJ=CMJ1*SIN (PH I RAD) 

CL  =  0 . 0 
GO  TO  900 


C  C 

C  INITIALIZE  TABLE  OF  EXPERIMENTAL  AND  CALCULATED  DATA  FOR  USE  C 

C  IN  INTERPOLATION  OF  FORCE  AND  MOMENT  VALUES  OF  CRUCIFORM  C 

C  CONFIGURATION  C 

C  C 

200  IF  (OPFLAG.NE .999)  CALL  TABULA  (OPFLAG) 

XJO  =  ( X JI N-XLE ) /LF 

IF  (XJO. GT.O. 19. AND.XJO.LT. 0.89)  GO  TO  230 
IF  (XJO. GE. 0.89)  GO  TO  220 

C  C 

C  IF  INPUT  VALUE  OF  XJ  IS  LESS  THAN  THE  MINIMUM  TABLE  VALUE*  C 

C  CALCULATE  FORCE  AND  MOMENT  USING  THE  LOWER  BOUND  TABLE  ONLY  C 

C  C 

XLO=. 1 9 
XHI=.19 
I BOUND  =  1 


CALL  INT (PHIJIN, PHIMIN,  ALPHIN  * DCMX 1 »DCMX 1 »  XLO ♦ XHI *  I  BOUND*  XJO*  CM J) 
CALL  INT (PHI JIN, PHIMIN ♦ ALPHIN *DCNX1 *DCNX1 ♦ XLO* XHI* IBOUND* XJO* CN) 
CMJ=CMJ*CT 
CN=CN*CT 

I F ( JUMPMP.NE . 0 ) GOTO  350 

CALL  INT (PHIJINfPHlMIN* ALPHIN* DC YX1* DC YX1»XL0»XHI, I80UND* X JO* CY ) 
CALL  INT (PHI JIN, PHIMIN, ALPHIN, DCHX1,DCHX1, XLO, XHI *I80UND, XJO, CNJ) 
CY=CY*CT 
CNJ=CNJ*CT 

CALL  INT (PHI JIN, PHIMIN, ALPHIN, DCLX1,DCLX1, XLO, XHI, IBOUND, XJO, CL) 


CL=CL*CT 
GO  TO  900 

C  C 
C  IF  INPUT  VALUE  OF  XJ  IS  GREATER  THAN  THE  MAXIMUM  TABLE  VALUE,  C 
C  CALCULATE  FORCE  AND  MOMENT  USING  THE  UPPER  BOUND  TABLE  ONLY  C 
C  C 


17 


220  XL0=.89 
XHI=.89 
IBOUND  =  1 

CALL  I NT (PHI JIN*PHIMIN*ALPHIN*DCMX3*DCMX3*XL0»XHI * I BOUND * XJO* CM J) 
CALL  INT <PHIJIN,PHlMIN»ALPHIN*DCNX3»DCNX3,XLO,XHI* IBOUND* X JO* CN) 
CMJ=CMJ*CT 
CN=CN*CT 

IF ( JUMPMP .NE.0) GOTO  350 

CALL  INT (PHI JIN,PHIM IN* ALPHIN*DCYX3*DCYX3*XL0*XHI* IBOUND* X JO* CY) 
CALL  INT (PHIJIN,PHIMIN* ALPH IN* DCHX3* DCHX3* XLO ♦ XH I ♦ I BOUND* XJO, CNJ) 
CY=CY*CT 
CNJ=CNJ*CT 

CALL  INT(PHIJIN,PHIMIN,ALPHIN.DCLX3*DCLX3,XL0*XHI*IB0UND*XJ0*CL) 


CL=CL*CT 
GO  TO  900 

C  C 

C  IF  INPUT  XJ  IS  WITHIN  BOUNDS*  AND  BETWEEN  MIDDLE  AND  UPPER  C 

C  VALUES.  INTERPOLATE  USING  MIDDLE  AND  UPPER  TABLES  C 

C  C 

230  XLO  =  .54 
XHI  =  .89 

IF  (XJO. GT. 0.19. AND. XJO.LE. 0.54)  GO  TO  240 


CALL  INT (PHIJIN.PHIMIN* ALPH IN* DCMX2* DCMX3* XLO* XHI , I BOUND » X JO* CM J) 
CALL  INT (PHI JIN,PHIMIN*ALPHIN,DCNX2»DCNX3*XL0*XHI . I  BOUND* XJO* CN) 
CMJ=CMJ*CT 
CN=CN*CT 

IF ( JUMPMP . NE . 0 ) GOTO  350 

CALL  INT (PHIJIN.PHIMIN* ALPHIN,DCYX2*DCYX3*XLO*XHI, IBOUND* XJO* CY) 
CALL  INT ( PHI  JIN, PHIM IN *ALPHIN*DCHX2*DCHX3.XLO* XHI* IBOUND* XJO. CNJ) 
CY=CY*CT 
CNJ=CNJ*CT 

CALL  INT (PHI JIN, PHIM IN. ALPH IN, DCLX2*DCLX3*XL0,XHI .IBOUND* XJO* CL) 


CL=CL*CT 
GO  TO  900 

C  C 

C  IF  INPUT  XJ  IS  WITHIN  BOUNDS,  AND  BETWEEN  LOWER  AND  MIDDLE  C 

C  VALUES*  INTERPOLATE  USING  LOWER  AND  MIDDLE  TABLES  C 

C  C 

240  XLO  =  .19 
XHI  =  .54 


CALL  INT (PHI JIN.PhIMIN, ALPHIN.DCMXl *DCMX2*XL0,XHI *  I  BOUND, XJO* CM J) 
CALL  INT (PHI JIN, PHIMIN. ALPH IN, DCNX1*DCNX2*XL0, XHI, IBOUND* X JO, CN) 
CMJ=CMJ*CT 
CN=CN*CT 

IF  < JUMPMP. NE. 0) GOTO  350 

CALL  INT ( PHIJIN.PHIMIN, ALPHIN*DCYX1 *0CYX2,XL0,XHI, IBOUND* XJO* CY) 
CALL  INT (PH I JIN, PHIM IN* ALPH IN*DCHX1»DCHX2» XLO, XHI, IBOUND* XJO* CNJ) 
CY=CY#CT 
CNJ=CNJ*CT  • 

CALL  INT (PH I JIN, PHIMIN ♦ ALPH IN*DCLX1 »DCLX2*XL0*XHI , IBOUND* XJO* CL) 


CL=CL*CT 
GO  TO  900 

C  C 
C  FOR  MONOPLANE  CONFIGURATIONS*  CHOOSE  THE  QUADRANT  IN  WHICH  C 
C  INPUT  QUANTITY  PHI  JET  APPEARS  C 
C  C 


oooooo  oooo  oooo  non  noo  ooo  ooo 


300  IF  (PHIJIN.GE. 0.0. AND. PHIJIN.LE. 90.0)  GO  TO  330 
IF  (PHIJIN. GT. 90.0. AND. PHIJIN.LE. 180.0)  GO  TO  320 
IF  (PHIJIN. GT. 180.0. AND. PHIJIN.LE. 270.0)  GO  TO  310 

IF  PHI  JET  IS  IN  QUADRANT  FOURt  ADJUST  ITS  INPUT  VALUE 

PHI JIN=PHI JIN-360 . 

JUMPMP  =  4 
GO  TO  340 

IF  PHI  JET  IS  IN  QUADRANT  THREE*  ADJUST  INPUT  VALUES  PHIJET  PHIM 

310  IF (PHIMIN.LT . 180. ) PH IM I N= 180. -PHIM IN 
IF (PHIMIN.GT.180. >PHIMIN=540.-PHIMIN 
PHIJIN=180. -PHIJIN 
JUMPMP=3 
GO  TO  340 

IF  PHI  JET  IS  IN  QUADRANT  TWO*  ADJUST  INPUT  VALUES  PHIJET*  PHIM 

320  IF (PHIMIN.LT. 180.) PHIM IN=PHIMIN* 100. 
IF(PHIMIN.GE.180.)PHIMIN=PHIMIN-180. 

PHIJIN=PHIJIN-180. 

JUMPMP=? 

60  TO  340 

IF  PHI  JET  IS  IN  QUADRANT  ONE*  ADJUST  INPUT  VALUES  PHIJET*  PHIM 

330  PHIMIN=360.-PHIMIN 
PHI JIN=-PHI JIN 
JUMPMPs 1 

INTERPOLATE  ADJUSTED  PHI  JET,  PHI  M  AND  ALPHA  VALUES 
NEEDED  FOR  EACH  QUADRANT 

340  GO  TO  200 

INITIALIZE  VARIABLES  USED  IN  THE  BODY-ALONE  PORTION  OF 
CALCULATIONS  FOR  THE  MONOPLANE  CONFIGURATION 

350  GO  TO  100 

355  IF  ( JUMPMP-2 )  360*370,380 

CALCULATE  FINAL  VALUFS  FOR  MONOPLANE  CONFIGURATION 
BY  ADJUSTED  QUADRANT  LOCATION  OF  PHIJET 

QUADRANT  ONE 

360  CY=-CNT*SIN(PHIRAD)»CT 
CNJ=-CMJ1*S IN (PH I RAD) 

CL=-.6*CN*SIN(PHIRAD) 

GO  TO  900 
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QUADRANT  TWO 

370  CMJ=-CMJ 
CN=-CN 

CY=-CNT*S IN (PHIRAD) *CT 
CNJ=-CMJ1*SIN< PHIRAD) 

CL=-0.6  *CN*SIN (PHIRAD) 

GO  TO  900 

QUADRANT  THREE 

380  IF  < JUMPMP.EQ.4)  GO  TO  390 
CM J=-CM J 
CN=-CN 

C Y=CNT*S IN (PHIRAD) *CT 
CNU=  CM J1*SIN (PHIRAD) 

CL= . 6*CN*SIN ( PHIRAD ) 

GO  TO  900 

QUADRANT  FOUR 

390  CY=CNT*SIN (PHIRAD) *CT 
CNJ=CMJ1*SIN (PHIRAD) 

CL=(0.6)*CN*SIN(PHIRAD) 

GO  TO  900 

WRITE  UNRECOVERABLE  ERROR  MESSAGE 
400  WRITE (5»410)NCONFG 

410  FORMAT ( 1X*33HC0NFIGURATI0N  NUMBER  IS  INCORRECT*  15/) 
900  RETURN 
END 
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SUBROUTINE  TA8ULA  (OPFLAG) 


THIS  SUBROUTINE  IS  CALLED  FROM  SUBROUTINE  SIDJET  TO  INITIALIZE 
DATA  MATRICES  USED  IN  INTERPOLATION  OF  NORMAL  Y  AND  Z  FORCES 
AND  MOMENTS 


4 

* 

4 

« 

4 

4 

» 

* 

4 

4 

* 

4 

4 

4 

ft 

4 

4 

4 

4 

* 

4 

4 

4 

4 

4 

4 

4 


COMMON  /MASTBL/  DCNX1 (13*3*7) .DCNX2 ( 13*3*7) .DCNX3 ( 1 3» 3« 7 > ♦ 
DCYX1 (13.3«7)  .DCYX2 (13.3,7) »DCYX3< 13.3,7) ♦ 
DCMX1 (13.3.7) *DCMX2( 13.3.7) .DCMX3 ( 1 3. 3. 7) , 
DCHX1 ( 13.3.7) *DCHX2( 13.3.7) .DCHX3< 13.3.7) . 
DCLX I (13.3.7) .DCLX2 (13.3.7) .DCLX3 ( 13*3.7) 
DIMENSION  AN  113(7) *AYI13(7) .AM 113(7). AH 113(7) »AL 113(7) ♦ 
AN213 (7) .AY213 (7) .AM213I7) .AH213(7) .AL213(7) . 

AN3 13 (7) .AY313 (7) .AM313(7) .AH313(7) *AL313(7)  . 

AN 123(7) *AY123(7) . AMI 23 (7) .AH123(7) .AL123(7) . 
AN223(7> * AY223 ( 7 ) .AM223(7) .AH223(7) *AL223(7)  . 
AN323 (7) ♦ AY323 (7) »AM323(7) .AH323(7) .AL323(7)  * 

AN 133(7) *AY133(7) .AM 133(7) .AH 133(7) .AL133(7) ♦ 
AN233 ( 7 ) ♦ AY233 ( 7> .AM233(7) .AH233(7) .AL233(7)  . 
AN333 ( 7 ) ♦ AY333 ( 7 ) .AM333(7) .AH333(7) .AL333(7> ♦ 

AN 153 (7) .AY153(7) *AM153(7) .AH 153 (7) ,AL153(7> ♦ 
AN253 (7) . AY253 ( 7) . AM253 ( 7 ) .AH253(7) .AL253(7) . 
AN353 ( 7 ) .AY353 (7) .AM353(7) ,AH353<7> .AL353(7) . 

AN  163 (7) .AY163(7) .AMI  63 (7) .AH163(7) .AL163(7) , 
AN263 ( 7 ) ♦ A Y263 ( 7 ) .AM263(7) .AH263(7) .AL263(7) ♦ 
AN363 ( 7 ) .AY  363 ( 7 ) .AM363(7) .AH363(7) *AL363(7) 
DIMENSION  AN  112(7). AY112(7). AM  112(7) .AH112(7).AL 112(7). 

AN312 (7) .AY312 (7) *AM312(7) .AH312(7) .AL312(7) . 
AN122 (7) .AY  122 (7) *AM122(7) .AH122(7) .AL122(7) . 
AN322 ( 7 ) ♦ AY322 (7) .AM322(7) »AH322(7) .AL322(7)  . 

AN  132(7) .AY132(7) .AM 132 (7) *AH132(7) »AL132(7) ♦ 
AN332 ( 7 ) ♦ A Y332 ( 7 ) .AM332(7) .AH332(7) .AL332(7) . 
AN152(7) .AY152(7) . AM 152 ( 7 ) ♦ AH 152 (7) *AL152(7) . 
AN352 ( 7 ) ♦ AY352 ( 7) .AM352(7) *AH352<7) .AL352 (7)  . 

AN 162(7) .AY  162 ( 7 ) . AMI  62 ( 7 ) . AH  162 (7) .AL162(7) ♦ 
AN362 ( 7 ) ♦ AY362 ( 7) .AM362(7) ,AH362(7) *AL362(7) 
INTEGER  OPFLAG 


INSERT  KNOWN  PARAMETER  VALUES  FROM  TABLES 

DATA  AN1 13/-. 0826.. 06.. 10*. 15*. 2*. 24. -.03/ 

DATA  AY11 3/. 41. .38 ..37*. 36. .40. .41.. 4 2/ 

DATA  AMI 13/-. 5.-1 . 9. - 1 . 9. - 1 . 93.-2. 2. -2. 2.- 1 . 4/ 

DATA  AH1 13/ 1.5* 1.6. 1.6. 1.66. 1.66. 1.8 .2. 6/ 

DATA  AL1 13/.0  25.  •  136*.  065*0.01  .-.02»”-.  125  .-.43/ 

DATA  AN213/ .2. .17.. 16*. 15. .15. .12*. 05/ 

DATA  AY213/.55».54*.52*.52..55..59».67/ 

DATA  AM213/-2.63*-2.63.-2.63.-2.63.-2.63»-2.63.-2.00/ 
DATA  AH21 3/3. 0.2. 8. 2. 6. 2. 57, 2. 3. 2. 3. 2. 6/ 

DATA  AL21 3/. 41*. 24 ..12. 0.08. 0.04 ,-.09. -.38/ 

DATA  AN313/-.85»-.8,-.8»-.8,-.75.-.7.-.5/ 

DATA  A Y31 3/ 1.1 6, 1.16. 1.16 ,1.1 6. 1.16. 1.16. 1.16/ 

DATA  AM313/-1. 8.-1. 8,-1. 8,-1. 8.-1. 8, -1.8. -1.8/ 

DATA  AH313/.24,.24..24,.24..24».24, .24/ 

DATA  AL31 3/. 00,. 080 ,.02,. 02*. 05.. 06*. 28/ 

DATA  AN1 23/. 13,. 25,. 15,. 15.. 2.. 2,. 02/ 

DATA  AY12 3/. 26. .4. .44. .36. .36, .39, .47/ 


21 


o  o  o  ooono 


DATA  AM123/-1 .3,-2.3,-1.93»-1.93»-2.2,-2.2.-l .6/ 
DATA  AH  123/2. 6, 1.6, 1.5,1.66,1.9.2.2,2.5/ 

DATA  AL123/.2, .15,.09,0.01 »-.08»-.16,-.4l/ 

DATA  AN223/.32».25».15».15».l7».13»0.0/ 

DATA  AY223/.3, .45* .50, .52, .58, .62, .68/ 

DATA  AM223/-3.0,-2.8,-2.7*-2.63»-2.7»-2.7,-2.5/ 
DATA  AH223/ 3. 0.2. 6 .2. 3 .2. 57 .2. 3. 2. 2 .2. 5/ 

DATA  AL223/. 30. .25.. 12. 0.08. 0.02. -.05. -.28/ 

DATA  AN 323/ -.7 .-.7 ♦ -.79*-. 8.” .73. -.73.- .55/ 

DATA  A Y323/ 1.1. 1.15 .1.155. 1.16. 1.16. 1.15* 1.1 4/ 
DATA  AM323/-2.0.-1 . 8.-1 .8 . -1 .8 .“1 . 8.-2. 0.-2. 0/ 
DATA  AH 323/. 2. . 15. .1». 24 ..28. .32. .32/ 

DATA  AL323/. 07 ,.07,0.02, 0.02, 0.02,0.02*. 24/ 

DATA  AN133/.35,.33,.24».15».17».12,-.07/ 

OATA  AY133/0.0..31..39* .36 ♦ . 33 ♦ . 34 * .48/ 

DATA  AM133/-2.0.-2.2,-2.0,”1.93,-1.8,-l.6,-.9/ 
DATA  AH  133/2. 9. 1.8* 1.3. 1.66.2.0.2.2.2.8/ 

DATA  AL 133/. 325 *.125. .06, 0.01. -.10*-. 19.-. 47/ 
DATA  AN2 33/. 42.. 26, .20.. 15, .15. .14. .13/ 

DATA  AY233/»20..40,.46..52».58..65».84/ 

DATA  AM233/-3.2,-2.9»-2.7.-2.63*-2.8.-3.0»-2.8/ 
DATA  AH 233/3. 2. 2.6.2. 57.2.57. 2. 4.2.2, 2.0/ 

DATA  AL233/. 24. .17. .15,0. 08, 0.01,-. 01.”. 03/ 

DATA  AN333/-. 66 • 76. - • 78 ♦ 8. -• 76 »- • 72. - • 56/ 
DATA  A Y333/ 1.14, 1.16, 1.16 ,1.16, 1.16. 1.16, 1.1 4/ 
DATA  AM333/-1.8.-1 .8.- 1 .8, *1 .8.-1 .8.-1 .8.-1 .8/ 
DATA  AH333/.24, .24. .24. .24. .24 ..24. .24/ 

DATA  AL333/.07. .040. .020,0.02.0.0,0.0» .23/ 

DATA  AN 153/. 4. .25. .2, . 15, .03, -.08 .-.44/ 

DATA  AY153/0.0..22».29».36..26».16,.43/ 

DATA  AM153/-2.4»-2.5,-2.5,-1.93.-1.6»-1.0,-.5/ 
DATA  AH 153/3. 1.2. 2. 1.66, 1.66.2.2.2.4.2.1/ 

DATA  AL 153/. 065. .026, .014,0.01 ,-.06,-. 14.-. 16/ 
DATA  AN253/. 22 ..18,. 16,. 15.. 05 .0.0, -.06/ 

DATA  AY253/.21. .40. .45, .52, .50, .60. .96/ 

DATA  AM253/-2.9»-?.8,-2.7,-2.63,-2.6,-2.5,-2.4/ 
DATA  AH 25 3/2. 8, 2. 6,2. 57,2. 57 »2. 6.2. 2. 1.6/ 

DATA  AL253/-. 15. 0.04 .0.08, 0.08.. 05 .0.08.. 3/ 

DATA  AN353/-.8,-»8»-»8.-»8.”»8.-.8.”.8/ 

OATA  AY353/ 1.1 2, 1.14, 1.15, 1.16, 1.06, 1.08, 1.1 3/ 

DATA  AM353/-1 .8.-1 .8, - 1 .8 ♦- 1 . 8 »— 1 . 8 » - 1 .8 » - 1 . 8/ 
DATA  AH353/.24. .24, .24. .24. .24, .24. .24/ 

DATA  AL353/0. 0,0. 0,0. 02. 0.02, 0.0 .0.02. .16/ 

DATA  AN163/.47, .4,. 28. . 15. .06, -.06,-. 38/ 

DATA  AY  163/0.0. .38, .4, .36, .29 ,.22. .28/ 

DATA  AM163/-3.7 .-3.2 »-2.8 .-1.93.- 1.93, -1.6. -1.0/ 
DATA  AH 163/4. 7, 2. 4 ,2.0 ,1.66, 2. 4, 2. 8. 3. 2/ 

DATA  AL 163/-. 4, 0.01,0.0 1,0. 01,0.0, -.08 ,.09/ 

DATA  AN263/.10. .18..20. .15. .08 ,0.0, -.23/ 

DATA  AY 263/ .52. .52.. 52. .52.. 48. .54, .72/ 

DATA  AM263/-2.9,-3.0,-2.8,“2.63,-2.8»-2.7,-2.2/ 
DATA  AH 263/ 2. 8,2. 7, 2. 6,2. 57, 2. 8, 2. 57.2.2/ 

DATA  AL263/-.3. .02. . 08. 0. 08, . 085. .17. .39/ 

DATA  AN363/-.6»-.7,-.7»-.8*-.8»-.8»-.95/ 

DATA  AY363/ 1.1 7. 1.19, 1.17. 1.16, 1.12, 1.08 ,1.05/ 
DATA  AM363/-1 . 8 ♦ - 1 .8 ♦ - 1 • 8 , ” 1 • 8 ♦- 1 .8. - 1 .8.-1 • 8/ 
DATA  AH363/.24».24..24..24..24..24»*24/ 
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DATA  AL363/.10* .10t . 08*0. 02. 0.0 *0. 0*0.0/ 

DATA  AN1 12/. 5 ,0.0 *-.15, -.15,-. 08, -.02* -.25/ 

DATA  AY112/-.1*.05,.10».15».17,.20,.29/ 

DATA  AMI 12/-4. 0,-2. 6, -2. 0,-2. 0,-2. 4, -2. 6, -2. 5/ 

DATA  AH1 12/2. 5, 2. 0,2. 0,2. 0,2. 0,2. 0*2. 5/ 

DATA  AL112/.15».04,.01»0.0*.06,.ll»-.23/ 

DATA  AN312/-1. 14,-1. 14,-1. 14,-1. 14,-1. 14,-1. 14 ,-1.14/ 
DATA  AY312/ 1.14* 1,14* 1.14* 1.14* 1.14* 1,14* 1.14/ 

DATA  AM312/-1.0»-1.0,-1.0»-1.0»-1.0*-1.0»-1.0/ 

DATA  AH31 2/ 1.0* 1.0 *1.0 *1.0* 1.0 *1.0 *1.0/ 

DATA  AL31 2/. 14* .055*. 0275* 0.0*-. 02* -.03* -.08/ 

DATA  AN122/.5,.02,-.15*-.l5,0.0»-.05*-.25/ 

DATA  AY  122/-. 2*. 07 *.15 *.15*. 15*. 21 *.38/ 

DATA  AM122/-3.6*-?.7»-2.0*-2.0*-2.2*”2.2*-2.0/ 

DATA  AH  122/2.8*1.6*1. 6*2. 0*2.0*2. 1*2.5/ 

DATA  AL12 2/ -.01*. 06*. 03*0.0* .035*. 07*-. 22/ 

DATA  AN322/-1. 14 *-l.l4»-l. 14,-1. 14,-1. 14,-1. 14 *-1.14/ 
DATA  AY322/ 1.14* 1.14* 1.14* 1.14*1.14*1.14*1.14/ 

DATA  AM322/-1.0*-1.0*-1.0*-1.0*-1.0*-1.0*-1.0/ 

DATA  AH322/ 1.0*1. 0*1. 0*1. 0*1. 0*1. 0*1.0/ 

DATA  AL322/. 04* .015*. 08*0.0* 0.0*0. 0*0.0/ 

DATA  AN 132/. 32* -.02* -.13* -.15 *-.08 *-.14 *-.3/ 

DATA  AY  132/-. 34, 0.0*. 12*. 15*. 06*. 11 *.26/ 

DATA  AM132/-3.0*-2.6*-2.0*-2.0*-2.0*-2.2*-2.2/ 

DATA  AH1 32/2. 6 *2. 2 *1.7 *2. 0*2. 2*2. 2 *2. 2/ 

DATA  AL1 32/ 0.0, 0.0, 0.0, 0.0,-. 01 »-. 02 *-.06/ 

DATA  AN332/-1.14»-1.14*-1.14*-1.14»-1.14*-1.14*-1.14/ 
DATA  A Y332/ 1.14, 1.14 *1.14* 1.1 4* 1.1 4, 1.14 *1.1 4/ 

DATA  AM332/ -1.0, -1. 0,-1. 0*-l»0*-1.0*-l. 0,-1. 0/ 

DATA  AH332/1, 0*1. 0*1. 0*1. 0*1. 0*1. 0*1.0/ 

DATA  AL332/ 0.0* 0.0 *0.0 *0.0 *0.0 *0.0 *0.0/ 

DATA  AN 152/0.0* -.10 *-.12.-. 15.-. 18. -.21 *-.3/ 

DATA  AY 152/-. 63*-. 04* . 15* . 15, .04 *-.02* .2/ 

DATA  AM152/-2.0*-2.0»-2.0»-2. 0,-2. 0,-2. 0* -2.0/ 

DATA  AH15 2/ 4. 0*2. 7*2. 0*2. 0*2. 0*2. 4*2.2/ 

DATA  AL 152/-. 17, -.07 *-.04.0.0*-. 06,0.0* . 17/ 

DATA  AN352/-1 .14.-1. 14,-1. 14,-1. 14 *-1.14*-l.l4*-l. 14/ 
DATA  A Y352/ 1.14, 1.14, 1.14, 1.14, 1.1 4. 1.14 ,1.1 4/ 

DATA  AM352/-1. 0,-1. 0,-1. 0,-1. 0,-1. 0,-1. 0,-1. 0/ 

DATA  AH352/1 .0, 1 .0, 1 .0, 1.0, 1.0, 1.0, 1.0/ 

DATA  AL352/-. 18, -.07 *-.04 *0.0,. 015*. 03, -.05/ 

DATA  AN162/. 1,0.0. -.05, -.15, -.12, -.08*. 04/ 

DATA  AY 162/-. 10* . 10*. 15* .15, .06* 0.0*-. 06/ 

DATA  AM162/-4.0*-3.0*-2.6*-2.0,-2.6»-2.6,-4.0/ 

DATA  AH162/4.4,2.7*2.2*2.0»2.6*3.0*4.0/ 

DATA  AL 162/-. 15* .08* .05*0.0, .05* .08* .2/ 

DATA  AN362/- 1.14 *-1.14, -1.1 4,-1. 14, -1.14,- 1.14,- 1.14/ 
DATA  AY362/ 1.14, 1.14, 1.14* 1.1 4, 1.1 4, 1.14, 1.14/ 

DATA  AM362/-1. 0,-1. 0,-1. 0,-1. 0,-1. 0,-1. 0,-1. 0/ 

DATA  AH362/ 1.0, 1.0, 1.0 *1.0, 1.0, 1.0, 1.0/ 

DATA  AL362/-. 15, -.04, -.02* 0.0, 0.0, .02*. 15/ 

EQUATE  THESE  TABLES  WITH  THEIR  POSITION  IN 
data  MATRICES 

DO  1  L  =  1 » 7 

DCNX1 (1,3,L)=AN113(L) 
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1  CONTINUE 
C 
C 

DO  2  L=  1  ♦  7 


2  CONTINUE 
C 
C 

DO  3  L=1  »7 


3  CONTINUE 
C 


OCYX1 ( 1 ♦ 3»  L ) =AY 1 1 3 (L ) 
DCMXl (I,3»L)=AM113(L) 
DCHX1 (1*3»L)=AH113(L) 
DCLX 1 ( 1  *  3»L ) =AL 1 1 3 (L ) 
DCNX2(1*3»L)=AN213(L) 
DCYX2(1.3»L)=AY213(L) 
DCMX2(lt3»L>=AM213(L) 
DCHX2(1»3»L)=AH213(L) 
DCLX 2 ( 1»3»L) =AL213(L) 
DCNX3(1»3»L)=AN313(L) 
DCYX3(1»3,L>=AY313(L) 
DCMX3(1»3»L)=AM313(L) 
DCHX3 ( 1  *  3*L) =AH31 3 (L ) 
DCLX3<1»3»L>=AL313(L> 


DCNX1 (2»3*L)=AN123(L) 
DCYX1 (2»3»L)=AY123(L) 
DCMXl (2t3»L)=AMl23(L) 
DCHX 1 ( 2»  3»  L ) =AH1 23 (L ) 
DCLX1 (2»3»L) =AL123(L) 
DCNX2(2»3»L)=AN223(L) 
DCYX2(2*3»L) =AY223(L) 
DCMX2(2»3»L)=AM223(L) 
DCHX2(2»3»L)=AH223(L) 
DCLX2(2»3»L)=AL223(L) 
DCNX3(2»3»L) =AN323(L) 
DCYX3(2»3»L)=AY323(L> 
DCMX3(2»3*U=AM323(L) 
DCHX3(2»3»L)=AH323(L) 
DCLX3(2»3»L)=AL323(L) 


DCNX1 (3*3»L)=AN133(L) 
DCYX1 (3*3»L) =AY133(L> 
DCMXl (3*3*L)=AM133<L) 
DCHXl (3»3»L)=AH133(D 
DCLX 1 ( 3»3»L ) =AL 1 33 (L ) 
DCNX2 (3»3*L)=AN233(L) 
DCYX2(3»3»L)=AY233(L) 
DCMX2(3»3.L)=AM233(L) 
DCHX2(3*3»L)=AH233(L) 
DCLX 2 ( 3»  3»L  >  =AL233 (L ) 
DCNX3(3»3»L)=AN333(L) 
DCYX3(3*3»L)=AY333(L) 
DCMX3(3»3»L)=AM333(L) 
DCHX3(3»3»L)=AH333(L) 
DCLX3(3»3*L)=AL333(L) 


C 

C 


C 
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c 

DO  4  L=1 *7 

DCNX1(5»3*L)=AN153(L) 
DCYX1 (5»3»L)=AY153(L) 
DCMX1 (5»3»L)=AM153(L) 
DCHX1 (5»3*L)=AH153(L> 
DCLX1 (5»3»L)=AL153(L) 
DCNX2(5.3»L)=AN253<L) 
DCYX2(5»3»L)=AY253(L) 
DCMX2<5.3*L)=AM253(L> 
DCHX2(5»3»L)=AH253(L) 
DCLX2(5,3»L)=AL253(L) 
DCNX3(5,3»L)=AN353«L) 
DCYX3(5,3»L)=AY353(L) 
DCMX3(5.3»L)=AM353(L) 
DCHX3(5»3*L)=AH353(L) 
DCLX3 (5»  3»L ) =AL353 (L ) 

4  CONTINUE 
C 

C 

DO  <5  L  =  1  ♦  7 

DCNX1 (6.3*L)=AN163(L) 
DCYXI (6»3»L)=AY163(L) 
DCMX1 (6«3»L>=AM163(L) 
DCHX1 (6*3»L)=AH163(L) 
DCLX1 (6»3»L)=AL163<L) 
DCNX2(6.3»L)=AN263(L) 
DCYX2(6»3*L)=AY263(L> 
DCMX2 ( 6* 3»L) =AM263  <L ) 
DCHX2(6»3»L)=AH263 (L) 
DCLX2(6»3»U  =AL263(L) 
DCNX3(6»3»L>=AN363(L) 
DCYX3(6»3»L) =AY363(L> 
DCMX3(6»3»L)=AM363(L) 
DCHX3(6*3»L)=AH363(L) 
DCLX3(6,3»L)=AL363(L) 

5  CONTINUE 


C 


C 

C 


DO  6  L  =  1 » 7 

DCNX1 (1*2*L)= AN 112(E) 
DCYX1 (1,2»L)=AY112(L) 
DCMX1  (1,2»L)=AMH2(L) 
DCHX1(1»2»L)= AH 112(1) 
DCLX1 (lt2»L)=AL112(L) 
DCNX3(1»2»L)=AN312(L) 
DCYX3(1*2*L)=AY312(L) 
DCMX3(1»2»L)=AM312(L) 
DCHX3(1»2»L)=AH312(L) 
DCLX3(1»2*L)=AL312(L) 

6  CONTINUE 
C 


C 
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DO  7  L  =  1  *  7 

DCNX 1 ( 2*2»L ) =AN122 (L ) 
DCYX1 (2*2»L) =AY122 (L) 
DCMX1 (2.2,L)=AM122(L) 
DCHX1 (2»2»L>=AH122(L) 
DCLX1 (2»2«L)=AL122<L) 
DCNX 3 (2»2»L) = AN 322 (L) 
DCYX3(2»2»L) =AY322(L> 
DCMX3(2*2«L)=AM322(L) 
DCHX3(2*2*L) =AH322(L) 
DCLX3(2»2»L)=AL322(L) 

7  CONTINUE 
C 

C 

DO  8  L= 1  *  7 

DCNX1 (3*2»L)=AN132(D 
DCYX1 (3»2tL)=AY132(L) 
DCMX1 (3»2»L)=AM132(L) 
DCHX1 (3»2»L)=AH132(L) 
DCLXI (3»2»L)=AL132(L) 
DCNX3(3«2»L)=AN332(L) 
DCYX3(3,2»L)=AY332(U 
DCMX3(3*2»L)=AM332(L) 
DCHX3(3.2»L)=AH332(L) 
DCLX3 ( 3»  2  *L ) =AL332 (L ) 

8  CONTINUE 
C 

C 

DO  9  L= 1 ♦ 7 

DCNX1 (5»2»L)=AN152(L) 
DCYX1 (5»2»L)=AY152(L) 
DCMXI (5«2»L)=AM152(D 
DCHX1 (5»2»L)=AH152(L) 
DCLX1 (5»2»L>  =AL152(L) 
DCNX3(5»2»L)=AN352(L) 
DCYX3(5.2»L)=AY352(L> 
DCMX3(5«2»L)=AM352(L) 
DCHX3 ( 5»  2»L  >  =  AH352 (L) 
DCLX3(5»2»L)=AL352(L) 

9  CONTINUE 
C 

C 

DO  10  L  =  1 » 7 

DCNX1 (6,2»L)=AN162(L) 
DCYX1 (6»2»L)=AY162(L) 
DCMXI (6»2.L) =AM162(D 
DCHX1 (6t2»L)=AH162(L) 
DCLXl (6*2»L) =AL162(L) 
0CNX3(6,2*L)=AN362<L> 
DCYX3(6»2»L)=AY362(L) 
DCMX3(6.2*L)=AM362(L) 
DCHX3(6»2»L) =AH362(L) 
DCLX3(6*2»L) =AL362 (L ) 

10  CONTINUE 
C 


C 

C 
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USE  THE  INITIAL  TABLE  ENTRIES  TO  CALCULATE  REMAINING 
VALUES  TO  FILL  DATA  MATRIX 

DO  12  J=l,3 
DO  11  L  =  1 » 7 

DCNX2(J»2*L) =DCNX1 <J»2»L> 
DCYX2(J»2,L)=DCYX1 (J.2,L) 

DCMX2( J.2»L) =DCMX1 (J»2*L) 
DCHX2(J»2*L)=DCHX1 ( J»  2»L ) 

DCLX2 ( J»2»L) =DCLX1 (J*2,L) 

11  CONTINUE 

12  CONTINUE 

DO  14  J=5»  6 
DO  13  L  =  l»7 

DCNX2 (J»2»L)=DCNX1 (J»2*L) 
DCYX2(J»2»L)=DCYX1 (J.2»L) 

DCMX2 ( J»  2»L ) =DCMX1 (J,2,L) 
DCHX2(Jt2*L)=DCHXl (J*2,L) 
DCLX2(J,2,L)=DCLX1 (J»2,L) 

13  CONTINUE 

14  CONTINUE 

DO  15  L  =  1 » 7 

DCNX1 (5»1*L)=-DCYX1 (1,3»L) 

DCYX1  ( 5» 1 *L ) =-DCNX 1 ( 1 ♦ 3*L) 
DCMX1(5*1»L)=-DCHX1 (1,3*L) 

DCHX1 <5* 1 *L) =-DCMX 1 ( 1  *  3»L) 

DCLX1 (5* 1 >L) =~DCLX1 ( 1  *  3  *  L ) 
DCNX2(5«1*L)=-DCYX2(1»3»L) 
DCYX2(5»1*L)=-DCNX2<1.3»L) 

DCMX2(5»1»L) =-DCHX2(l»3»L) 
DCHX2<5»1»L)=-DCMX2(1,3,L) 
DCLX2(5»1»L)=-DCLX2(1,3»L) 
DCNX3(5»1»L)=-DCYX3(1»3»L) 
DCYX3(5tl*L)=-DCNX3(l«3*L) 
DCMX3(5»1*D=-DCHX3(1»3*L) 

DCHX3 (5» 1 »L ) =-DCMX3 ( 1 » 3»L) 
DCLX3(5»1,L)=-DCLX3(1»3*L) 

15  CONTINUE 

DO  16  L=1 ♦ 7 

DCNX1 (7*3»L) =OCNXl (1»3»8-L) 

DCYX1 (7*3.L)=DCYX1 (1,3,8-L) 

DCMX 1 ( 7  *  3  *  L ) =DCMX 1 ( 1 ♦ 3»  8-L ) 

DCHX1 (7»3»L)=DCHX1(1»3,8-L) 

DCLX1 ( 7»  3*L ) =DCLX 1 ( 1 ♦ 3*8-L ) 

DCNX2 ( 7*  3»L ) =DCNX2 ( 1 ♦ 3»8-L) 
DCYX2(7*3*L>=DCYX2(1*3»8-L) 

DCMX2(7*3»L) =DCMX2(1»3.8-L) 

DCHX2(7*3»L) =DCHX2(1.3»8-L) 

DCLX2(7»3»L) =DCLX2 ( 1 *3»  8-L ) 

DCNX3 ( 7*  3»L) =DCNX3 ( 1 »3,8-L) 
DCYX3(7»3*L)=DCYX3(1»3»8-L> 

DCMX3(7»3»L) =DCMX3(1*3»8-L) 

DCHX3(7*3»L) =DCHX3 ( 1 ♦ 3»  8-L) 

DCLX3(7»3»L) =DCLX3(l»3t8-L) 


27 


non 


16  CONTINUE 

DO  17  L=1 *7 

DCNX1 < 11*1*L) =-DCYXl <1*3*8-L) 
OCYX1 ( lltl»L) =-DCNXl (1*3* 8-L) 
OCMX1 ( 1 1*1»L)=-DCHX1 <1*3*8-L) 
DCHX 1 ( 1 1 » 1 * L ) =-DCMX 1 ( 1 ♦ 3 » 8-L ) 
DCLX1 ( 11»1.L) =-DCLXl <1»3*8-L) 
DCNX2 ( 1 1  *  1 »L ) =-DCYX2 ( 1  *  3*8-L ) 
DCYX2(11»1*L)=-DCNX2<1»3*8-L) 
DCMX2<11*1*L) =-DCHX2(l*3*8-L) 
DCHX2( 1 1*1*L) =-DCMX2(l»3*8-L) 
0CLX2<11*1*L)=-DCLX2<1»3*8-L) 
DCNX3(11*1»L) =-DCYX3 ( 1  *  3»  8-L ) 
DCYX3 ( 1 1  *  1 »L  > =-DCNX3 ( 1*3* 8-L ) 
DCMX3 ( 1 1 » 1 *L) =-0CHX3 ( 1 » 3*8-L) 
DCHX3 ( 1 1 ♦ 1 » L ) =-DCMX  3  <  1 * 3 ♦ 8-L ) 
DCLX  3 ( 1 1 ♦ 1 *L ) =-DCLX3 ( 1*3*  8-L ) 

17  CONTINUE 

DO  18  L  =  1 *  7 

0CNX1(4*1»L)=-DCYX1 (2»3*L) 
OCYX1 (4»1»L)=-0CNX1 (2*3*L) 
DCMX1 (4* 1*L) =-DCHXl (2*3*L) 
DCHX1 (4» 1»L) =-OCMXl (2t3*L) 
OCLX1 (4* 1*L) =-OCLXl <2*3*L> 
OCNX2 (4*1 *L)=-DCYX2(2*3*L) 
DCYX2(4* 1*L>=-DCNX2(2*3»L) 
0CMX2(4*ltL) =-DCHX2(2»3»L) 
DCHX2(4*1*L>=-DCMX2(2»3*L) 
DCLX2(4* 1»L)=-0CLX2(2»3»L) 
DCNX3<4*1*L>  =-DCYX3<2*3*L) 
DCYX3<4* 1*L> =-0CNX3(2t3»L) 
DCMX3(4*1*L)=-DCHX3(2»3*L) 
DCHX3(4»1*L) =-0CMX3<2»3»L) 
DCLX3(4»1»L) =-DCLX3 (2»3*L) 

18  CONTINUE 

DO  19  L  =  1 *  7 

DCNX1 <8*3*L) =DCNX1 (2*3*8-L) 
OCYX1 <8*3»L) =DCYX1 <2* 3. 8-L) 
OCMX1 (fl*3tL) =DCMX1 (2* 3* 8-L) 
OCHX1 (8»3»L)=0CHX1 (2»3»8-L) 
DCLX 1 ( 8  *  3*L ) =DCLX 1 (2* 3* 8-L) 
DCNX2<8*3*L) =DCNX2(2*3«8-L) 
DCYX2(8»3*L)=DCYX2(2*3*8-L) 
DCMX2 (8*3*L) =DCMX2 (2*3»8-L> 
DCHX 2 ( 8  *  3*L ) =DCHX2 ( 2*3*  8-L) 
DCLX2  <8*3*L) =DCLX2(2*3*8-L> 
DCNX3 ( 8*  3*L  > =DCNX3 ( 2*3*  8-L ) 
DCYX3(8*3*L)=DCYX3(2»3*8-L) 
DCMX3 (8*3»L) -DCMX3 (2*3»8-L) 
DCHX3(8*3»L>  =DCHX3(2*3*8-L) 
DCLX3(8*3*L) =DCLX3 ( 2*3*  8-L ) 

19  CONTINUE 

DO  20  L  =  1 » 7 

DCNX1 (10* 1»L) =-OCYXl (2*3»8-L) 
DCYX1 ( 10* 1*L) =-DCNXl (2, 3* 8-L) 
DCMX1 ( 10* 1*L) =-DCHXl (2, 3* 8-L) 
DCHX 1 ( 1 0  *  1 *L) =-DCMX 1 (2, 3* 8-L) 


DCIX1 ( 1 0* 1 *L) =-DCLXl (2,3*8-L) 
OCNX2(10» 1*L) =-DCYX2(2,3»8-L> 
DCYX2(10.1»U=-DCNX2(2,3»8-L) 
DCMX2 ( 10*1 ♦  !_) =-DCHX2 (2»3»8“L) 
DCHX2(10*1»L) =-DCMX2 (2»  3*  8-1) 
DCLX2 ( 1 0  *  1 *L) =-DCLX2 (2*  3»8-L) 
DCNX3(10»1»L) =-DCYX3(2*3»8-L> 
DCYX3(10* ltL) =-DCNX3(2,3»8-L) 
DCMX3(10» 1,L) =-DCHX3(2.3»8-L) 
DCHX3(10»1»L) S-DCMX3 ( 2*  3*8-L) 
0CLX3 ( 1 0 1 1 »L) =-DCLX3 ( 2»  3*  8“L) 

20  CONTINUE 

DO  21  L=1.7 

DCNX1(3«1»L)=-DCYX1 <3*3tL> 
DCYX1(3»1»L) =-DCNX I ( 3»3*L ) 
DCMX1(3»1»L)=-DCHX1 (3*3»L) 
DCHX1(3»1*L)=-DCMX1 (3*3»L> 
0CLX1<3*1*L)=-DCLX1 <3.3*L) 
DCNX2(3* 1»L)=-DCYX2(3»3»L) 
DCYX2(3» 1*L)=-DCNX2(3»3*L) 
DCHX2(3* l*L)=-DCHX2(3»3tL) 
DCHX2(3»1 »L) =-DCMX2(3»3*L) 
DCLX2(3*1»L)=“DCLX2(3*3»L) 
DCNX3(3« 1*L)=-DCYX3(3»3»L) 
DCYX3  <3*1 »L)=-DCNX3(3»3»L> 
DCMX3(3» 1*L)=-DCHX3<3»3*L) 
DCHX3(3»1»U=-DCMX3(3»3»L> 
DCLX3(3* 1*L) =-DCLX3(3»3»L) 

21  CONTINUE 

DO  22  L= 1 » 7 

DCNX1 (9»3tL)=DCNXl <3»3t8-L) 
DCYX1 (9*3»L) =DCYX1 (3»3»8-L) 
DCMX1 (9*3»L) =DCMX1 (3»3,8-L> 
DCHX1 (9»3*L)=DCHX1 <3»3*8-L) 
DCLX1 ( 9 « 3  »L ) =DCLX 1 <3*3,8-L> 
DCNX2(9»3»L) =DCNX2(3*3»8~L) 
DCYX2(9*3*L) =DCYX2(3*3»8-L) 
DCMX2(9*3»D=DCMX2(3*3»8-L) 
DCHX2(9*3»L>=DCHX2(3»3t8-L> 
DCLX2(9»3*L) =DCLX2(3»3»8-L) 
DCNX3(9,3*L)=DCNX3(3»3t8-L) 
DCYX3(9»3»L)=DCYX3(3»3»8~L) 
DCMX3(9»3»L) =DCMX3 (3,3,8-L) 
DCHX3(9»3»L) =DCHX3 ( 3* 3,8-L ) 
DCLX3(9»3»L)=0CLX3(3»3»8-L) 

22  CONTINUE 

DO  23  L=1 »7 

DCNX1 (9» 1»L) =-DCYXl (3*3»8-L) 
DCYX1 (9» 1»L)=-DCNX1 (3»3»8-L) 
DCMX1 (9*1 *L) =-DCHXl (3»3*8-L) 
DCHX1  (9*1*U=-DCMX1  <3*3»8-L> 
DCLX1 (9»1»L)=-DCLX1 (3»3»8-L) 
DCNX2(9,1,L)=-DCYX2(3.3»8-L) 
DCYX2 ( 9» 1 »L)=-DCNX2(3»3»8-L) 
DCMX2C9*1*L)*-DCHX2(3«3*8-L) 
DCHX2(9» 1»L)=-DCMX2(3»3»8-L) 
DCLX2(9*l»L)=-DCLX2(3t3*8-L) 
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UCNX3 I9t 1*L) =-DCYX3t3»3*8-L) 
DCYX3(9»1»L) =-DCNX3 ( 3*  3*  8-L ) 
DCMX3(9»1»L)=-DCHX3<3*3»8-L) 
DCHX3<9» 1*L) =-DCMX3(3»3*8-L) 
DCLX3(9»1»L) =  -DCLX3 ( 3*  3»  8-L ) 

23  CONTINUE 

DO  24  L=1 *  7 

DCNX1 (7»I»L)=-DCYX1 (5»3»8-L) 
DCYX1 ( 7» 1 »L ) =-DCNX 1 (5»  3 ♦ 8-L ) 
DCMX I ( 7 ♦ 1 *L ) =-DCHX 1 <5* 3* 8-L) 
DCHX1 ( 7* 1 *L ) =-DCMX 1 (5 1 3»8-L) 
DCLX1 (7»1»L)=-DCLX1 <5»3»8-L> 
DCNX2 ( 7» 1 »L ) =-DCYX2  <5*  3 1 8-L ) 
DCYX2 ( 7* 1 *L ) =-DCNX2 (5*3»  8-L ) 
DCMX 2 ( 7  *  1 *L ) =-DCHX2 (5*  3*  8-L ) 
DCHX2(7»1»L)=-DCMX2(5»3»8-L> 
DCLX2  < 7* 1 tL ) =-DCLX2 (5t  3»8-L ) 
DCNX3 ( 7  *  1 »L ) =-DCYX3 (5»  3*  8-L ) 
DCYX3(7»1»L) =-DCNX3 (5*  3* 8-L ) 
DCMX3 ( 7* 1 *L ) =-DCHX3 (5»3*8-L ) 
DCHX3(7» I*L) =-DCMX3(5»3*8-L) 
DCLX3(7*1*L) =-DCLX3  <5*  3*  8-L ) 

24  CONTINUE 

DO  25  L=  1  *  7 

DCNX 1 ( 1 1 ♦ 3 ♦ L  > =DCNX 1 ( 5 ♦ 3 ♦ 8-L ) 
DCYX 1 ( 1 1 ♦ 3*L ) =DCYX 1 (5.3»8-L> 
DCMX1 ( 1 1 *3*L) =DCMX1 <5»3»8-L> 
DCHX1 (11,3*L) =DCHXI (5«3»8-L) 
DCLX1 (11»3»L)=DCLX1 (5»3*8-L) 
DCNX2( 11»3»L) =DCNX2(5»3*8-L) 
DCYX 2 ( 1 1 *3*L) =DCYX2 (5*3»8-L) 
DCMX2(11*3*L)=DCMX2 <5»3*8-L) 
DCHX2( 1 1»3»L) =DCHX2(5»3*8-L) 
DCLX2(11»3»L)=DCLX2(5»3*8-L) 
DCNX3U  I.3»L)=DCNX3(5»3»8-L) 
DCYX3(11*3,L)=DCYX3(5*3»8-L) 
DCMX3( 11»3»L) =DCMX3(5»3»8-L) 
DCHX3 ( 1 1 1 3»L ) =DCHX3 (5*  3*  8-L ) 
DCLX3(11*3*L)=DCLX3(5*3*8-L) 

25  CONTINUE 

DO  28  L= 1 ♦ 7 

DCNX1 (It I»L)=-DCYX1 (5»3»L> 
DCYX1 ( 1»1»L) =-DCNXl <5»3»L) 
DCMX1 ( 1»1»L) =-DCHXl (5»3»L> 
DCHX 1 ( 1  *  1 »L) =-DCMXl (5*3»L) 
DCLX 1 ( 1 1 1 *L  > =-DCLX 1 (5»  3»L) 
DCNX2(1 »1»L> =-DCYX2(5»3*L> 
DCYX 2  < 1  *  1 »L) =-DCNX2 (5*3»L) 
DCMX 2 ( 1  *  1 tL ) =-DCHX2 (5»  3*L) 
DCHX2(1.1*L)=-DCMX2(5*3*L) 
DCLX2 ( 1 ♦ 1 *L) =-DCLX2  <5t  3»L) 
DCNX 3 ( 1  *  1  *  L  > =-DCYX3 (5«  3»L) 
DCYX3 ( 1 ♦ 1 *L) =-DCNX3 (5»3»L) 
DCMX3 ( 1  *  1 *L) =-DCHX3 (5*  3»L) 
DCHX3( 1*1»L> =-DCMX3(5»3*L) 
DCLX3 ( 1 ♦ 1 *L) =-DCLX3 (5»  3» L) 
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28  CONTINUE 

DO  27  K=1 *3 

DO  26  L= 1 ♦ 7 

DCNX1 (13*K,L)=DCNX1 ( 1 »K«L> 
DCYX1 (13,K,L) =DCYX1 ( I*K*L) 
DCMX] < 13tK*L) =DCMX1 ( 1 *K*L) 
DCHX1 (13»K*L) =DCHX 1 ( I *K  *L ) 
OCLX1 (13»K,L)=DCLX1 <1»K*L) 
DCNX2(13*K,L)=DCNX2(1,K»L) 
DCYX2(13»K,U  =DCYX2(1*K»L> 
DCMX2( 13*K,L) =DCMX2(1.K»L) 
DCHX2(13»K,L)=DCHX2(1*K»L> 
DCLX2<13.K,L)=DCLX2U,K»L) 
DCNX3(13»K»L) =DCNX3 ( I*K*L) 
DCYX3<13»K»L)=DCYX3(1»K»L) 
DCMX3(13*K,L) =DCMX3(1»K,L) 
DCHX3<  13*  K»L)=DCHX3(1  ♦«*»_) 
DCLX3(13,K,L) =DCLX3(1»K»L) 

26  CONTINUE 

27  CONTINUE 

DO  29  L=1 *7 

DCNX1(6»1»L)=~DCYX1 (6*3»8-L) 
DCYX1 (6»1»L)=-DCNX1 (6»3»8-L) 
DCMXl(6fl*L)=-DCHXl (6.3i8-L) 
DCHX1(6»I*L)=-DCMX1 <6*3»8-L> 
DCLX1(6.1»L)=-DCLX1 <6»3*8-L) 
DCNX2(6»1*L)=-DCYX2(6.3*8-L) 
DCYX2(A,UL)=-DCNX2(6»3»8-L) 
DCMX2(6» 1*L)=-DCHX2(6,3»8-L) 
DCHX2 ( 6 » 1 *L ) =-DCMX2(6.3»8-L) 
DCLX2(6«1*L) =-DCLX2 (6*3*8-L) 
DCNX3(A*1»L)=-DCYX3(6*3*8-L) 
DCYX3(6*l»L>=-DCNX3(6»3t8-L) 
DCMX3(6» 1 tL)=-DCHX3(6*3»8-L> 
DCHX3 ( 6 ♦ 1 *  L ) =-DCMX3 ( 6  »  3*  8-L ) 
DCLX3(6* 1»L)=-DCLX3(6»3»8-L) 

29  CONTINUE 

DO  30  L  =  1 » 7 

DCNX 1 ( 12*  3*L) =DCNXl (6»3»8~L) 
DCYX1 (12»3»L)=DCYX1 (6»3»8-L> 
DCMX1 (12*3*L)=DCMX1 <6.3»8-L> 
DCHX1 (12»3.L) =DCHX1 (6*3.8-L) 
DCLX1 (]2.3»L) =DCLX1 (6»3»8-L) 
DCNX2<I2*3»L) =DCNX2(6»3.8-L) 
DCYX2<12»3»L)=DCYX2(6*3»8~L) 
DCMX2(12»3.L) =DCMX2(6»3*8-L) 
DCHX2(I2«3»L)=DCHX2(6»3*8"L) 
DCLX2 ( 1 2 ♦ 3*L) =DCLX2 (6 » 3»8”L ) 
0CNX3<12»3»L)=DCNX3(6»3»8-L> 
DCYX3(12*3*L) =DCYX3(6»3*8-L) 
DCMX3< 12»3»L) =DCMX3(6»3»8-L) 
DCHX3< 12*3tL) =DCHX3 <6.3*8«L> 
DCLX3(12»3.L)=DCLX3(6»3*8-L) 


30  CONTINUE 

DO  31  L=1 »7 

OCNX1 <12»1»L)=-DCYX1 <6»3»L> 
DCYX1 ( 12»1*L) =-OCNXl <6*3»L> 
0CMX1(12»1,L)=-DCHX1 (6»3»L) 
OCHX1 ( 12. l.L) =-OCMXl (6»3»L) 
DCLX1 (12,1,L)=-DCLX1(6,3»L) 
DCNX2(]2*1,L>  =-OCYX2(6«3*L) 
DCYX2( 12»1,L)=-DCNX2(6»3»L) 
DCMX2(12*1»L) =-0CHX2 (6» 3*L ) 
0CHX2(12.1*L)=-0CMX2(6,3.L) 
DCLX2(12»1,L)=-0CLX2(6»3»L) 
DCNX3( 12*1 ,L)=-DCYX3 <6,3, L) 
DCYX3( 12»1*L) =-DCNX3(6»3»L) 
OCMX3( 12»1»L) =-DCHX3(6,3»L> 
0CHX3( 12.1»L)=-0CMX3(6.3»L) 
DCLX3 ( ) 2* 1 ♦ L) =“DCLX3 (6»  3»L) 

31  CONTINUE 

DO  32  L=l»7 

OCNX1 ( 7 » 2*L ) =DCNX1 ( 1 *2t 8-L) 
DCYX1 (7*2*L)=DCYX1 ( 1 *2*8-L) 
DCMX1 (7»2*D  =DCMX 1 ( 1 t 2*8-L) 
DCHX1 (7»2*L) =DCHX1 ( 1»2,8-L) 
DCLX1 (7»2»L)=DCLX1 (1*2»8-L) 
0CNX2 ( 7 ♦ 2*L ) =DCNX2 ( 1 »2*8-L) 
DCYX2(7«2»L) =DCYX2(1»2«8-L) 
DCMX2(7»2»L) =DCMX2 ( 1 ♦2»8-L) 
DCHX2(7»2*L) =DCHX2(1.2»8-L) 
0CLX2(7.2»L) =DCLX2(1»2»8-L) 
DCNX3(7*2*L) =DCNX3( 1»2»8-L) 
DCYX3(7*2*L) =DCYX3( 1*2«8-L) 
DCMX3  <  7  *  2*L ) =DCMX3 ( 1 *2»8-L ) 
0CHX3(7»2,L) =DCHX3( 1*2*8-L) 
0CLX3(7»2*L)=0CLX3(1»2»8-L) 

32  CONTINUE 

00  33  L= 1  *  7 

DCNX1 ( 11,2,L) =-OCYXl (1*2*8-L) 
DCYX1 ( 1 1  *  2»L ) =-OCNXl ( 1 «2*8-L ) 
DCMX1 < 1 1 *2»L) =-DCHXl (1»2»8-L) 
0CHX1 (11»2,L)=-DCMX1(1»2*8-L) 
DCLX1 ( 1 1 *2»L) =-OCLXl ( 1 *2»8-L) 
DCNX2(11*2»L) =-0CYX2 ( 1 ♦ 2»8-L ) 
DCYX?(11»2*L) =-0CNX2 ( 1 »2*8"L ) 
DCMX2 ( 1 1 »2» L) =-DCHX2 (1*2*  8-L ) 
DCHX2( 1 1*2*L) =-DCMX2<l*2*8-L> 
DCLX2<11*2*L)=-DCLX2<1*2*8-L> 
DCNX3<11*2*L)=-DCYX3<1*2*8-L> 
DCYX3< 11*2«L)=-DCNX3(1*2*8-L) 
DCMX3<11*2*L)=-DCHX3<1*2*8-L> 
DCHX3 ( 1 1 *2»L) =-DCMX3< 1*2*8-L> 
DCLX3 ( 1 1 »2*L> =-DCLX3 ( 1 *2* 8-L) 

33  CONTINUE 

DO  34  L  =  1 » 7 

DCNX1 <4*2*1 >=-DCYXl (2*2 *l> 
DCYX1 ( 4  *  2  *  L ) =-OCNXl <2.2»L> 
DCMX1 <4*2*L) =-DCHXl <2»2*L> 
DCHX1 (4*2*L)=-DCMX1 (2*2*L» 
DCLX1 (4*2*1) =-DCLXl <2*2*L) 


DCNX2(4t2»L> =-DCYX2 < 2  *  2* L ) 
DCYX2(4»2*L)=-DCNX2(2»2»L) 
0CMX2 ( 4 ♦ 2»L) =-DCHX2 <  2*  2*L ) 
DCHX2(4.2*L)=-0CMX2(2»2»L> 
DCLX2(4»2*L)=-DCLX2(2*2*L) 
DCNX3(4»2*t)=-DCYX3(2*2*L> 
DCYX3(4»2»L>  =-DCNX3(2.2*L) 
0CMX3(4»2»L)=-DCHX3(2»2»L) 
DCHX3(4»2»L)=-DCMX3(2»2»L) 
0CLX3(4«2»L)=-DCLX3<2»2*L) 

34  CONTIMUE 

DO  35  L=1  »7 

OCNX1(0»2,L)=DCNX1 (2.2»8-L) 
DCYX1 (8»2.L)=DCYX1 (2,2,8-L) 
DCMX1 (8*2»L)=0CMX1 (2»2»8-L) 
DCHX1 (8»2»L)=DCHX1 (2»2»8-L) 
DCLX1 (8«2*L)=DCLX1(2*2»8-L> 
DCNX2(8*2»L) =0CNX2(2»2»8-L) 
0CYX2(8»2*L)=DCYX2(2»2»8-L) 
DCMX2<8»2»L)=DCMX2(2.2,8-L) 
DCHX2(8*2»L)=DCHX2(2»2,8-L) 
0CLX2(8»2*L)=DCLX2(2,2»8-L) 
DCNX3(8t2»L)=DCNX3(2*2*8-L) 
DCYX3(8i2»L>=DCYX3(2»2»8-L> 
DCMX3 ( 8»2»L)=DCMX3(2*2«8-L) 
DCHX3(8»2»L>=0CHX3(2*2*8-L) 
DCLX3(8»2,L) =DCLX3 < 2 ♦ 2, 8-L ) 

35  CONTINUE 

00  36  L  =  1 » 7 

DCNX1 ( 10*2fL) =-DCYXl (2»2»8-L) 
OCYX 1 ( 1 0 ♦ 2»L ) =-DCNX 1 (2»2»8-L) 
DCMX 1 ( 1 0  *  2  *L) =-DCHXl 1 2»  2  *  8~L) 
DCHX1 (}0»2»L)=-DCMX1 (2»2»8-L) 
DCLX1 (]0»2»L)=-DCLX1 <2»2»8-L> 
OCNX2(10»2*L)=-DCYX2(2»2*8-L) 
0CYX2( 10*2.L) =-DCNX2(2,2«8-L) 
DCMX2<10»2»L) =-DCHX2 ( 2*2*  8-L) 
DCHX2(10.2fL)=-OCMX2(2»2«8-L) 
0CLX2( 10»2»L) =-DCLX2(2*2*8-L> 
OCNX.3(10*2»L)=-OCYX3(2»2»8-L) 
OCYX3(10»2»L)=“OCNX3(2»2*  8“L) 
DCMX3(10*2»L)=-DCHX3(2»2#8-L) 
DCHX3(10»2*L)=-DCMX3(2»2*8-L) 
DCLX3(10*2.L)=-DCLX3(2»2*8-L) 

36  CONTINUE 

DO  37  L= 1 »7 

DCNX1 <9i2*L>=DCNX1 (3»2*8-L) 
DCYX1 (9»2»L) =0CYX1 (3»2»8-L) 
DCMX1 <9.2*L)=DCMX1 (3,2,8-L) 
OCHX 1 (9*2*L)=DCHX1 <3»2»8-L) 
DCIX1  <9»2*U=DCLX1  (3.2*8-L) 
DCNX2(9»2»L)=DCNX2(3»2»8-L) 
DCYX2(9»2*L)=DCYX2(3»2,8-L) 
0CMX2(9»2.L)=DCMX2(3,2t8>L) 
DCHX2(9»2»L) =0CHX2 ( 3»2*8-L ) 
DCLX2(9»2.L>=DCLX2(3,2»8-L> 
DCNX3(9»2*L)=DCNX3(3»2*8-L) 
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DCYX3(9*2»L>  =DCYX3(3*2*8-L) 
DCMX3(9»2.L>=DCMX3(3*2.8-L> 

DCHX3(9*2*L) =DCHX3(3»2*8“L) 
DCLX3(9*2.L)=DCLX3(3*2*8-L> 

37  CONTINUE 

DO  30  L=1  *7 

DCNX1 ( 12*2*L) =DCNXl (6.2*8-L> 

DCYX1 (l2f2»L)=DCYXl (6*2* 8-L) 

DCMX1 <12»2,L)=DCMX1 <6»2*0-L) 

DCHX1 (12*2.L)=DCHXl <6*2*8-L> 

DCLX1 (12*2,L)=DCLX1 <6»2*8-L) 

DCNX2<12*2*L>  =DCNX2<6.2.8-L) 
DCYX2(12*2*L)=DCYX2 (6*2*8-L ) 
DCMX2n2*2,L)=DCMX2<6*2*8-L) 
DCHX2(12.2*L)=DCHX2(6*2*8-L) 

DCLX2(I2,2,L) =DCLX2 (6,2.8-L) 
DCNX3(12*2*L)=DCNX3<6*2*8-L) 

DCYX3 ( 12»2*L) =DCYX3 (6*2*8-L) 
DCMX3(12*2*L)=DCMX3(6*2»8-L> 
DCHX3(12*2*L>=DCHX3<6»2*8-L) 

DCLX3 ( 1 2. 2.L) =DCLX3 <6* 2* 8-L) 

38  CONTINUE 

DO  39  L  =  I  ♦  7 

DCNX1 (2*1,L>=0.5*(DCNX1 <1*1»L> *DCNX1  <3*1*L) > 
DCYX1(2»1*L>=0.5*(DCYX1(1*1*L)*DCYX1  (3.1*L) > 
DCMX1 (2*1*L)=Q.5*(DCMX1 (1»1.L) ♦DCMXl (3.1*L) > 
DCHX1 (2*1»L)=0.5*(DCHX1 ( 1  *  1 *L) *DCHX1 <3,1*L>  > 
DCLXl <2.1 *L> =0.5* (DCLXl (1.1. L) ♦DCLXl (3.1.L) > 
DCNX2(2.1.L)=0.5*<DCNX2<1.1.L>  *DCNX2 <3. 1 .L)  ) 
DCYX2(2»1*L)=0.5*(DCYX2(1*1»L) *DCYX2 ( 3, 1 .L ) ) 
DCMX2 (2. 1  »L) =0 .5* (DCMX2  ( 1 » 1 .L) *DCMX2 ( 3. 1 *L ) ) 
DCHX2(2»1.L)=0.5*(DCHX2(1»1»L) ♦DCHX2 ( 3. 1 »L ) ) 
DCLX2(2,1,L>=0.5*<DCLX2<1,1,L>  *DCLX2 < 3. 1 .L ) ) 
DCNX3<2.1.L>=0.5*<DCNX3<1.1»L)*DCNX3<3.1.L) ) 
DCYX3(2.1»L>=0.5*(DCYX3<1.1.L> *DCYX3 ( 3, 1 *L ) ) 
DCMX3<2.1.L)=0.5*<DCMX3<1.1.L>*DCMX3<3.1.L>  ) 
DCHX3(2»1*L>=0.5*(DCHX3<1.1.L> *DCHX3 ( 3. 1 .L ) ) 
DCLX3  <2. 1  .L) =0 .5* (DCLX3 ( 1 ♦ 1 *L) ♦DCLX3 ( 3. 1 *L ) > 

39  CONTINUE 

DO  40  L=1 »7 

DCNX1 (4.3*L)=0.5*(DCNX1 <3»3.L)*DCNX1 (5.3.L) > 
DCYX1 <4.3.L)=0.5*(DCYX1 (3.3.L) +DCYX1 (5.3.L) > 
DCMX1 <4»3.L>=0.5*<DCMX1 (3.3.L) *DCMX1 (5,3.L) ) 
DCHX1 (4*3*L)=0.5*(DCHX1 (3*3*L) ♦DCHXl (5*3*L) ) 
DCLXl (4*3»L) =0.5* (DCLXl (3»3»L) *DCLX1 <5.3*L> ) 
DCNX2(4»3*L)=0.5*(DCNX2(3*3*L) +DCNX2 (5*3*L ) ) 
DCYX2(4»3*L)=0.5*(DCYX2(3,3*L) ♦DCYX2 (5*3*L) ) 
DCMX2(4*3*L)=0.5*(DCMX2(3*3*L) *DCMX2 (5*3*L) ) 
DCHX2(4*3*L)=0.5*(DCHX2(3*3»L) *DCHX2 < 5* 3»L ) > 
DCLX2(4*3*L)=0.5*(DCLX2(3»3*L) *DCLX2 (5.3.L) ) 
0CNX3(4»3*L)=0.5*(DCNX3(3»3*L) *DCNX3 (5.3.L) ) 
DCYX3(4.3.L)=0.5*<DCYX3<3.3.L> ♦DCYX3(5»3.L) ) 
DCMX3(4»3*L)=0.5*(DCMX3(3*3*L) *DCMX3 (5.3.L) ) 
DCHX3(4*3»L)=0.5*(DCHX3(3*3»L) *DCHX3 <5. 3.L > ) 
DCLX3(4«3*L>=0.5* (DCLX3 (3*  3*L) *DCLX3 (5*  3»L ) ) 
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40  CONTINUE 

DO  41  L=l*7 

DCNX1  (8t  1  *0=0*5*  (DCNX1  (7»  1  *L)  ♦DCNXl  (9*  1»L)  ) 
DCYX1 (8»1*L> =0.5* (DCYXl (7* 1 »L> ♦DCYXl (9»1»L> > 
DCMX1 (8»1*L)=0.5*(DCMX1 (7» 1 ♦DCMX1 (9» 1 *L) ) 
DCHX1 (8*1*L)=0.5* (DCHX 1 (7* 1 »L) ♦DCHX1 (9» 1 »L) ) 
DCLX1 (8*1*L>=0.5* (DCLX1 ( 7* 1 »L) ♦DCLXl (9* 1 »L> ) 
DCNX2 (8*1»DS0.5*  (DCNX2 ( 7* 1 »L)  ♦DCNX2 (9« 1 »L)  ) 
DCYX2(8*1»L)=0.5* (DCYX2 (7*1*L)*DCYX2(9»1»L) ) 
DCMX2<8*1»L)=0.5*(DCMX2<7,1»L)*DCMX2(9,1»L>) 
DCHX2 (8» 1 *L) =0.5* (DCHX2 (7» 1 *L) *008X2(9* 1 »L) ) 
DCLX2(8tl»L>=0.5*(DCLX2<7»l»L>*DCLX2<9»l»L> ) 
DCNX3 (8*1*L)=0.5* (DCNX3 ( 7 1 1 *L) +DCNX3 (9» 1 *L ) ) 
DCYX3 (8»1*L)=0.5* (DCYX3 (7*1»L)*DCYX3(9*1*L) ) 
DCMX3 (8»1*L)=0.5* (DCMX3 (7  »1»L)*DCMX3(9»1»L) ) 
DCHX 3 (8*1*L)=0.5*(DCHX3(7*1*L)*DCHX3(9*1*L) ) 
DCLX3(8»1»L)=0.5*(DCLX3(7*1»L)»DCLX3(9*1*L) ) 

41  CONTINUE 

DO  42  L=1  »7 

DCNX1 < 10»3»L) =0.5* (DCNXl (9»3*L) ♦DCNXl (11*3»L) ) 
DCMX1 <10»3.L) =0.5* (DCHX l(9»3tL) ♦DCMXl < 1 1 »3»L) ) 
DCYXl <10»3*L) =0.5* (DCYXl <9» 3*L> *DCYX1 (11 »3*L> ) 
DCHX1 ( 10»3»L) =0.5* (DCHXl (9»3»L) ♦DCHXl ( 11.3.L) ) 
DCLX1 ( 10»3»L> =0 .5* (DCLX 1 <9*  3*L) ♦DCLXl ( 11?3*L) ) 
DCNX2 (10.3*L)=0.5*(DCNX2(9«3»L) *DCNX2 ( 1 1 . 3.L ) ) 
DCYX2 (10»3»L)=0.5*(OCYX2(9»3»L)*OCYX2(11*3»L)) 
DCHX2 ( 10i3»L> =0.5* (DCMX2(9*3»L) ♦DCMX2 ( 11* 3*L> ) 
DCHX2 ( 1 0 . 3.L ) =  0 .5* (DCHX2  <9t  3*L) ♦DCHX2 ( 1 1 »3*L ) ) 
DCLX2(10*3*L)=0.5*(DCLX2(9»3»L)*OCLX2<11.3fL) ) 
DCNX3 ( 10»3»L)=0.5*(DCNX3(9»3*L) *DCNX3 ( 1 1 *3*L ) ) 
DCYX3(10t3*L)=0.5*(DCYX3(9»3*L)*DCYX3(ll»3*L) ) 
DCMX3 (10*3*L)=0»5* (DCMX3 (9»  3»L) ♦DCMX3(11*3»L) ) 
DCHX 3 (10»3»L)=0«5*(DCHX3(9»3*L)*DCHX3(11»3»L) ) 
DCLX 3 (10»3»L>=0«5* (DCLX3 (9t  3»L) ♦DCLX3(11»3*L) ) 

42  CONTINUE 
0PFLAG=999 
RETURN 
END 


ooo  ooo  oooo  o  o  ooo  oooooooo 


SUBROUTINE  INT (PHIK , PHI J. ALPHA, ADC1 ,ADC2, 

*  XLOtXHI * IBOUNO* AXIN* AXOUT ) 

THIS  SUBROUTINE  SEARCHES  A  RANGE  TABLE  TO  CHOOSE 
BRACKETING  VALUES  FOR  INPUT  VARIABLES*  THEN  INTERPOLATES 
THESE  VALUES  LAST-TO-F IRST 


FIND  BRACKETING  INDICES  AND  VALUES  FOR  PHIM  PHIJET  AND  ALPHA 

DIMENSION  RTJ ( 13) *RTK(3) *RTL<7) *ADC1 ( 13.3,7) ♦ ADC2 ( 1 3, 3. 7> 

SET  UP  RANGE  TABLES  FOR  PHIJET.  PHIM  AND  ALPHA 

DATA  RTJ/0. 0,22.5.45.0 ♦67.5,90.0.135. 0.180. 0,202. 5, 

*  225. 0,247. 5,270. 0f315. 0,360.0/ 

DATA  RTK/-22.5.-45.0.-67.5/ 

DATA  RTL/- 10.0, -4.0, -2. 0,0.0 ,2.0. 4.0, 10.0/ 

LOCATE  UPPER  AND  LOWER  VALUES  OF  PHIM.  AND  THEIR  CORRESPONDING 
SUBSCRIPTS;  SAVE  these  AS  NEW  VARIABLES 

DO  1  J=1 ♦ 12 

IF  (PHIJ.GE.RTJ(J)  .AND.PHIJ.LE.RTJ(JM)  )  GO  TO  10 

1  CONTINUE 

SAVE  VALUES  AND  SUBSCRIPTS  FOR  INTERPOLATION  ROUTINES 

10  JLO= J 
JH I = J* 1 
VALJLO=RTJ(J> 

VAL JHI=RT  J ( J* 1 ) 

DO  THE  SAME  FOR  PHIJET  AND  ALPHA 
15  DO  2  K=1 ♦ 2 

IF  (PHIK.LE.RTK(K) . AND. PHIK .GE .RTK (K* 1 ) )  GO  TO  20 

2  CONTINUE 
GO  TO  60 

20  KLO=K ♦ 1 
KHI  =K 

VALKHI=RTK(K) 

VALKLO=RTK ( K ♦ 1 ) 

C  C 

25  00  3  L=1 ,6 

IF  ( ALPHA, GE. RTL (L) . AND • ALPHA. LE. RTL (L*l) )  GO  TO  30 

3  CONTINUE 
GO  TO  70 

30  LLO=L 
LHI=L*1 
VALLLO=RTL (L ) 

VALLHI=RTL (L*  1 ) 

GO  TO  100 

C  C 
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ooo  ooo  oooo  o  o  ooo  oooooooo 


nno  n  o  onooo 


C  IF  either  phik  or  alpha  are  out  of  bounds*  set  these  equal  to  the 

C  MINIMUM  OR  MAXIMUM  RANGE  VALUES  AND  CONTINUE 
C 

60  IF  (PHIK. GT. -22.5)  GO  TO  65 
PHlK=-67.5 
GO  TO  15 
65  PHlK=-22.5 
GO  TO  15 
C 

70  IF  (ALPHA. GT. 10)  GO  TO  75 
ALPHA=-10.0 
GO  TO  25 
75  ALPHA=10.0 
GO  TO  25 


DO  ACTUAL  INTERPOLATION  GIVEN  PREVIOUSLY  CALCULATED 
BRACKETING  VALUES 

100  CALL  TRITER(JLO*JHl*KLO»KHI*LLO*LHl* 

*  ALPHA, PHIK*PHIJ»AOUTl, 

*  VALJHI , VALKHI , VALLHI ,VALJLO, VALKLO, VALLLO* ADC  1 ) 
IF  ( IBOUND.EQ. 1 )  GO  TO  900 

CALL  TRITER(JLO,JHl,KLO,KHl,LLO,LHl, 

*  ALPHA, PHIK*PHIJ,A0UT2, 

*  VALJHI, VALKHI, VALLHI, VALJLO, VALKLO* VALLLO, ADC2) 

DO  MATRIX-TO-MATRIX  FINAL  INTERPOLATION 

AXOUT=( ( (AXIN-XLO)* ( A0UT2-A0UT1 ) )/(XHI-XLO> ) ♦AOUT1 
GO  TO  999 
900  AXOUT  =  AOUT 1 
999  RETURN 
END 
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ooo  o  o  ooooo  o  ooo 


ooo  ooo  ooo 


SUBROUTINE  TRUER  ( JLO*JHI  ,KLO»KHI  *LLO*LHI » 

*  ALPHA»PHIK*PHIJ»AOUT» 

*  VALJHI *VALKHI *VALLHI *VALJLO* VALKLO* VALLLO* ADC) 

C  C 

DIMENSION  ADC ( 13*3»7) ,AL < 13*3) » AK < 13) 

INTERPOLATE  IN  ALPHA 

DO  5  J=JLO*JHI 
DO  4  K=KHI*KLO 

AL(J»K) =( ( (ALPHA-VALLLO)# (ADC( J*K,LHI ) -ADC ( J*K*LLO) ) ) 

*  / (VALLH I -VALLLO) ) *ADC ( J*K*LLO) 

4  CONTINUE 

5  CONTINUE 

INTERPOLATE  IN  PHIJET 
200  DO  6  J=JL0* JHI 

AK ( J) = ( ( (PHIK-VALKLO)*(AL(J*KHl)-AL(J»KLO) ) ) 

*  /(VALKHI-VALKLO) ) ♦AL(J»KLO) 

6  CONTINUE 

INTERPOLATE  IN  PHIM  LAST 

400  AOUT=< ( (PHIJ-VALJL0)*(AK(JHI)-AK(JL0) ) ) 

*  / ( VALJHI -VALJLO) ) ♦AK ( JLO) 

RETURN 
END 
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